GEOMETRICAL  MODELS  OF  THE  PHASE  SPACE  STRUCTURES 
GOVERNING  REACTION  DYNAMICS 


By 

Holger  Waalkens 

and 

Stephen  Wiggins 


IMA  Preprint  Series  #  2269 

(August  2009) 


INSTITUTE  FOR  MATHEMATICS  AND  ITS  APPLICATIONS 

UNIVERSITY  OF  MINNESOTA 

400  Lind  Hall 
207  Church  Street  S.E. 

Minneapolis,  Minnesota  55455-0436 

Phone:  612-624-6066  Fax:  612-626-7370 
URL:  http://www.ima.umn.edu 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

AUG  2009  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2009  to  00-00-2009 

4.  TITLE  AND  SUBTITLE 

Geometrical  Models  of  the  Phase  Space  Structures  Governing  Reaction 
Dynamics 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Minnesota, Institute  for  Mathematics  and  Its 

Applications, Minneapolis, MN, 55455-0436 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION 

_ _ _  OF  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

45 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Geometrical  Models  of  the  Phase  Space  Structures  Governing 

Reaction  Dynamics 

Holger  Waalkens^’^  and  Stephen  Wiggins^ 

June  26,  2009 


^  School  of  Mathematics,  University  of  Bristol,  University  Walk,  Bristol  BS8  ITW,  UK 

^  Department  of  Mathematics,  University  of  Groningen,  Nijenborgh  9,  9747  AG  Groningen,  The  Netherlands 
E-mail:  h.waalkens@math.rug.nl,  s.wiggins@bristol.ac.uk 

Abstract 

Hamiltonian  dynamical  systems  possessing  equilibria  of  saddle  x  centre  x  •  •  •  x  centre  stability  type 
display  reaction-type  dynamics  for  energies  close  to  the  energy  of  such  equilibria;  entrance  and  exit  from 
certain  regions  of  the  phase  space  is  only  possible  via  narrow  bottlenecks  created  by  the  influence  of 
the  equilibrium  points.  In  this  paper  we  provide  a  thorough  pedagogical  description  of  the  phase  space 
structures  that  are  responsible  for  controlling  transport  in  these  problems.  Of  central  importance  is  the 
existence  of  a  Normally  Hyperbolic  Invariant  Manifold  (NHIM),  whose  stable  and  unstable  manifolds  have 
sufficient  dimensionality  to  act  as  separatrices,  partitioning  energy  surfaces  into  regions  of  qualitatively 
distinct  behaviour.  This  NHIM  forms  the  natural  (dynamical)  equator  of  a  (spherical)  dividing  surface 
which  locally  divides  an  energy  surface  into  two  components  (‘reactants’  and  ‘products’),  one  on  either 
side  of  the  bottleneck.  This  dividing  surface  has  all  the  desired  properties  sought  for  in  transition  state 
theory  where  reaction  rates  are  computed  from  the  flux  through  a  dividing  surface.  In  fact,  the  dividing 
surface  that  we  construct  is  crossed  exactly  once  by  reactive  trajectories,  and  not  crossed  by  nonreactive 
trajectories,  and  related  to  these  properties,  minimizes  the  flux  upon  variation  of  the  dividing  surface. 

We  discuss  three  presentations  of  the  energy  surface  and  the  phase  space  structures  contained  in  it  for 
2-degree-of- freedom  (DoF)  systems  in  the  threedimensional  space  R®,  and  two  schematic  models  which 
capture  many  of  the  essential  features  of  the  dynamics  for  n-DoF  systems.  In  addition,  we  elucidate  the 
structure  of  the  NHIM. 


Contents 


1  Introduction  2 

2  “React ion- Type”  Dynamics  and  Phase  Space  Transition  State  Theory  3 

2.1  The  normal  form .  4 

2.2  Explicit  definition  of  the  phase  space  structures  in  the  normal  form  coordinates .  6 

2.3  The  foliation  of  the  reaction  region  by  Lagrangian  submanifolds .  11 

2.4  Implications  for  the  original  system .  11 

3  Models  of  the  local  phase  space  structures  for  2-DoF  systems  12 

3.1  Model  I:  projection  of  2-spheres  to  their  equatorial  planes .  14 

3.2  Model  II:  stereographic  projection  of  2-spheres .  21 

3.3  Model  III:  McGehee  representation .  26 

3.4  Implications  for  Nonlinear  Hamiltonian  Vector  Fields .  27 

4  A  model  of  the  n-DoF  system  in  the  space  of  the  integrals  of  motion  34 


1 


5  Structure  of  the  NHIM  for  n-DoF  35 

5.1  The  NHIM  (n  =  3  DoF) .  35 

5.2  The  NHIM,  forn  >  3  DoF .  38 

6  Conclusions  39 

7  Acknowledgements  41 


1  Introduction 

Transition  state  theory  (TST)  is  one  of  the  most  important  theoretical  and  computational  approaches  to 
analyzing  chemical  reaction  dynamics,  both  from  a  qualitative  and  quantitative  point  of  view.  Transition 
state  theory  was  created  in  the  1930’s,  with  most  of  the  credit  being  given  to  Eyring,  Polanyi,  and  Wigner, 
who  are  referred  to  as  the  “founding  trinity  of  TST”  in  Miller’s  important  review  on  chemical  reaction  rates 
([Miller(1998)]).  Nevertheless,  important  contributions  were  also  made  by  Evans,  Farkas,  Szilard,  Horiuti, 

Pelzer,  and  Marcelin,  and  these  are  described  in  the  discussions  of  the  historical  development  of  the  subject 
given  in  [Laidler  &  King(1983),  Poliak  &  Talkner(2005)]. 

The  central  ideas  of  TST  are  fundamental  and  provide  a  very  natural  and  fruitful  approach  to  “transfor¬ 
mation”  in  many  areas  far  removed  from  chemistry.  For  example,  it  has  been  used  for  certain  studies  in  atomic 
physics  [Jaffe  et  al.(2000)],  studies  of  the  rearrangements  of  clusters  [Komatsuzaki  &  Berry(1999)],  studies  of 
transport  phenomena  in  solid  state  and  semi-conductor  physics  [Jacucci  et  al.(1984),  Eckhardt(1995)],  stud¬ 
ies  of  diffusion  dynamics  in  materials  [Voter  et  al.(2002)],  cosmology  [de  Oliveira  et  al.(2002)],  and  celestial 
mechanics  [Jaffe  et  al.(2002),  Waalkens  et  al. (2005b)]. 

The  fundamental  assumptions  of  TST  we  clearly  stated  by  Wigner,  who  begins  by  stating  that  he  considers 
chemical  reactions  in  a  setting  where  the  equilibrium  Maxwell-Boltzmann  velocity  and  energy  distributions 
are  maintained  (see  [Mahan(1974)]  for  a  detailed  discussion  of  this  point)  and  for  which  the  potential  energy 
surface  is  known  ([Garrett (2000)]).  He  then  gives  the  following  assumptions  from  which  he  derives  TST: 

1.  the  motion  of  the  nuclei  occurs  on  the  Born-Oppenheimer  potential  energy  surface  (“electronic  adia- 
baticity”  of  the  reaction) 

2.  classical  mechanics  adequately  describes  the  motion  of  the  nuclei 

3.  there  exists  a  hypersurface  in  phase  space  dividing  the  energy  surface  into  a  region  of  reactants  and  a 
region  of  products  having  the  property  that  all  trajectories  that  pass  from  reactants  to  products  must 
cross  this  dividing  surface  precisely  once. 

It  is  important  to  note  that  Wigner  clearly  developed  his  ideas  in  phase  space,  the  arena  for  dynamics.  It 
is  important  to  keep  this  in  mind  since  a  great  deal  of  later  developments  occur  in  configuration  space,  in 
which  certain  dynamical  properties  are  obscured.  Nevertheless,  the  typical  mental  picture  for  the  geometry 
associated  with  TST  is  built  around  a  saddle  point  of  a  potential  energy  surface  for  a  one  degree-of-freedom 
(DoF)  or  two  DoF  system,  i.e.  it  starts  with  a  configuration  space  notion  of  sufficiently  low  dimension  that 
standard  visualisation  intuition  applies. 

It  follows  from  Wigner’s  assumptions  that  fundamental  geometric  structure  of  TST  is  a  dividing  surface 
which  locally  divides  the  energy  surface  into  two  disjoint  components  (which  we  refer  to  as  the  “bottleneck 
property” )  and  which  is  free  of  local  “recrossing”  (otherwise,  the  flux  through  the  dividing  surface  would  be 
overestimated).  These  properties  are  crucial  for  reaction  rate  computations. 

The  problem  of  how  to  define  and  construct  a  dividing  surface  with  these  properties  was  solved  for 
two  degrees  of  freedom  (DoF)  in  the  70’s  and  80’s  by  Pechukas,  Poliak  and  others  [Child  &  Pollak(1980), 
Pechukas(1981),  Pechukas  &  McLafferty(1973),  Pechukas  &  Pollak(1978),  Pechukas  &  Pollak(1977),  Pechukas  &  Pollak(l 
Poliak  et  al.(1980),  Pollak(198I),  Poliak  &  Pechukas(1979),  Poliak  &  Child(1980)]  who  constructed  the  di¬ 
viding  surface  from  a  periodic  orbit,  the  so  called  periodic  dividing  surface  (PODS). 

The  generalization  to  systems  with  more  than  two  degrees  of  freedom  has  posed  severe  problems. 

Periodic  orbits  lack  sufficient  dimensionality  for  the  construction  of  dividing  surfaces  for  systems  with 
three  or  more  degrees  of  freedom.  Progress  in  constructing  a  dividing  surface  for  systems  with  three 
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or  more  degrees  of  freedom  was  made  with  the  introduction  of  a  fundamentally  new  object,  a  nor¬ 
mally  hyperbolic  invariant  manifold,  or  “NHIM”  ([Wiggins(1990),  Wiggins(1994)]),  that  takes  the  place 
of  the  periodic  orbit.  The  NHIM  is  the  main  building  block  (and  the  natural  higher  dimensional  gener¬ 
alization  of  the  periodic  orbit  used  by  Pechukas,  Poliak  and  others)  of  the  phase  space  transition  state 
theory  mentioned  above.  This  has  been  described  more  fully  in  recent  work  has  shown  that  the  fun¬ 
damental  assumptions  of  transition  state  theory  can  only  be  realized  in  phase  space  for  systems  with 
three  or  more  DoF.  This  work  is  described  in  a  series  of  papers  [Wiggins  et  al.(2001)],[Uzer  et  al.(2002)], 
[Waalkens  et  al. (2004a)],  [Waalkens  &  Wiggins (2004)],  [Waalkens  et  al. (2004b)], [Waalkens  et  al. (2005a)], 
[Waalkens  et  al.  (2005c)],  [Schubert  et  al.(2006)],  [Waalkens  et  al.(2008)]  where  the  fundamental  dynamical 
and  geometrical  framework  for  phase  space  TST  is  developed.  Morever,  it  is  shown  that  new  types  of  “phase 
space  structures”  are  responsible  for  governing  the  dynamics. 

Chief  among  these  structures  is  the  NHIM.  ^  This  invariant  manifold  is  normally  hyperbolic,  i.e.,  it  is  of 
saddle  type,  having  stable  and  unstable  manifolds  which  are  of  sufficient  dimensionality  to  act  as  separatrices, 
i.e.  impenetrable  barriers  to  the  dynamics  which  serve  to  partition  energy  surfaces  into  regions  of  qualitatively 
distinct  behaviour. 

In  addition,  the  NHIM  forms  the  natural  (dynamical)  equator  of  a  (spherical)  dividing  surface  (DS) 
which  locally  divides  the  energy  surface  into  two  components,  one  on  either  side  of  the  bottleneck.  The 
dividing  surface  is  crossed  (locally)  exactly  once  by  those  trajectories  which  pass  from  one  component  to  the 
other;  it  is  free  of  local  recrossing  by  trajectories  and  is  of  minimal  directional  flux,  and  is  thus  the  optimal 
dividing  surface  sought  for  in  variational  transition  state  theory  [Waalkens  &  Wiggins(2004)].  In  a  nutshell, 
the  dividing  surface  enables  a  natural  definition  of  the  phase  space  regions  corresponding  to  reactants  and 
products.  In  chemical  terms,  the  equator  of  the  dividing  surface,  the  NHIM,  forms  the  energy  surface  of  the 
activated  complex,  consisting  of  an  oscillating  supermolecule  poised  between  the  reactants  and  the  products. 

The  purpose  of  this  paper  is  to  develop  the  geometrical  intuition  associated  with  these  high  dimensional 
phase  space  structures  that  govern  reaction  dynamics  by  constructing,  step-by-step,  three  representations  of 
the  phase  space  geometry  for  n  =  2  DoF,  and  another  model  that  illustrates  various  aspects  of  the  phase 
space  structures  for  an  arbitrary  number  n  of  DoF.  Moreover,  we  discuss  the  internal  structure  of  the  NHIM. 
Before  proceeding  to  the  develpment  of  these  models  we  give  a  technical  summary  of  the  relevant  phase  space 
structures,  their  interperation  in  terms  of,  and  implications  for,  reaction  dynamics,  and  the  method  by  which 
they  are  realised  (the  Poincare-Birkhoff  normal  form  theory) .  All  of  the  material  described  here  can  be  found 
in  the  above  references. 


2  “Reaction- Type”  Dynamics  and  Phase  Space  Transition  State 
Theory 

We  consider  a  Hamiltonian  system  with  n  degrees  of  freedom,  phase  space  coordinates  (Q,  P)  €  x  and 
Hamiltonian  function  H .  We  assume  that  (Qo?  ^o)  is  equilibrium  point  of  Hamilton’s  equations  which  is 
of  saddle-centre-. .  .-centre  stability  type.  ^  By  adding  a  constant  energy  term  which  does  not  change  the 
dynamics  we  can  achieve  H{Qq^  Pq)  =  0,  and  for  simplicity  of  exposition,  we  can  assume  that  the  coordinates 
have  been  suitably  translated  so  that  the  relevant  equilibrium  point  (Qo?  ^o)  is  at  the  origin.  For  much  of  the 
discussion  below,  we  will  consider  iso-energetic  geometrical  structures  belonging  to  a  single  positive  energy 

^The  concept  of  an  “invariant  manifold”  is  fundamental  to  the  description  and  quantification  of  a  wide  variety  of  dynamical 
phenomena.  However,  the  use  of  the  term,  and  its  significance,  may  not  be  so  clear  to  someone  with  a  non-dynamical  systems 
oriented  background.  Therefore  we  give  a  brief  explanation  of  the  term,  from  which  will  follow  its  implications  for  trajectories  in 
phase  space.  For  our  purposes  we  can  think  of  the  term  “manifold”  as  being  synonomous  with  the  term  “surface”.  “Invariance” 
of  the  manifold  can  be  characterized  in  two  (equivalent)  ways.  One  is  that  trajectories  through  initial  conditions  on  the  manifold 
must  be  contained  in  the  manifold  for  their  entire  pasts  and  futures.  The  other  way  of  characterizing  an  invariant  manifold  is 
by  saying  that  the  vector  field  giving  rise  to  the  dynamics  is  tangent  to  the  manifold  at  every  point  of  the  manifold.  It  follows 
immediately  from  both  characterizations  of  invariance  that  any  trajectory  with  initial  condition  not  on  the  invariant  manifold 
can  never  “cross”  the  invariant  manifold.  In  this  way,  invariant  manifolds  of  codimension  1,  i.e.  of  dimension  one  less  than  the 
dimension  of  phase  space,  act  as  separatrices.  They  form  the  “barriers  to  transport”  or  “impenetrable  barriers  in  phase  space”. 

^We  will  define  this  more  precisely  shortly.  However,  briefly,  it  means  that  the  matrix  associated  with  the  linearization 
of  Hamilton’s  equations  about  this  equilibrium  point  has  two  real  eigenvalues  of  equal  magnitude,  with  one  positive  and  one 
negative,  and  n  —  1  purely  imaginary  complex  conjugate  pairs  of  eigenvalues.  We  will  assume  that  the  eigenvalues  satisfy  a 
nonresonance  condition  that  we  will  describe  more  fully  in  the  following. 
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surface  S(/i)  :=  H~^{h)  for  some  constant  h  >  0.  In  addition,  we  will  typically  restrict  our  attention  to 
a  certain  neighbourhood  C,  local  to  the  equilibrium  point.  We  will  defer  until  later  a  discussion  of  exactly 
how  this  region  is  chosen;  suffice  it  to  say  for  now  that  the  region  is  chosen  so  that  a  new  set  of  coordinates 
can  be  constructed  (the  normal  form  coordinates)  in  which  the  Hamiltonian  can  be  expressed  (the  normal 
form  Hamiltonian)  such  that  it  provides  an  integrable  nonlinear  approximation  to  the  dynamics  which  yields 
phase  space  structures  to  within  a  given  desired  accuracy. 

Locally,  the  (2n  —  l)-dimensional  energy  surface  S(/i)  has  the  structure  of  x  R  in  the  2n-dimensional 
phase  space.  The  energy  surface  I](/i)  is  split  locally  into  two  components,  “reactants”  and  “products”  ,  by 
a  (2n  —  2)-dimensional  “dividing  surface”  that  is  diffeomorphic  to  and  which  we  therefore  denote  by 

S^"“^(/i).  The  dividing  surface  that  we  construct  has  the  following  properties:- 

•  The  only  way  that  trajectories  can  evolve  from  reactants  to  products  (and  vice-versa),  without  leaving 
the  local  region  £,  is  through  §j"“^(ft,).  In  other  words,  initial  conditions  (ICs)  on  this  dividing  surface 
specify  all  reacting  trajectories. 

•  The  dividing  surface  that  we  construct  is  free  of  local  recrossings;  any  trajectory  which  crosses  it  must 
leave  the  neighbourhood  C  before  it  might  possibly  cross  again. 

•  The  dividing  surface  that  we  construct  minimizes  the  flux,  i.e.  the  directional  flux  through  the  dividing 
surface  will  increase  upon  a  generic  deformation  of  the  dividing  surface  (see  [Waalkens  &  Wiggins  (2004)] 
for  the  details). 

The  fundamental  phase  space  building  block  that  allows  the  construction  of  a  dividing  surface  with  these 
properties  is  a  particular  Normally  Hyperbolic  Invariant  Manifold  (NHIM)  which,  for  a  fixed  positive  energy 
h,  will  be  denoted  SNfflM(^)-  The  NHIM  is  diffeomorphic  to  and  forms  the  natural  dynamical  equator 

of  the  dividing  surface:  The  dividing  surface  is  split  by  this  equator  into  (2n  —  2)-dimensional  hemispheres, 
each  diffeomorphic  to  the  open  (2n  —  2)-ball,  We  will  denote  these  hemispheres  by  and 

B‘^~^(h)  and  call  them  the  “forward  reactive”  and  “backward  reactive”  hemispheres,  respectively.  B'^f^ih) 
is  crossed  by  trajectories  representing  “forward”  reactions  (from  reactants  to  products),  while  B‘^~^{h)  is 
crossed  by  trajectories  representing  “backward”  reactions  (from  products  to  reactants). 

The  (2n  —  3)-dimensional  NHIM  can  be  viewed  as  the  energy  surface  of  an  (unstable)  invariant  subsystem 
which  as  mentioned  above,  in  chemistry  terminology,  corresponds  to  the  “activated  complex” ,  which  as  an 
oscillating  “supermolecule”  is  located  between  reactants  and  products. 

The  NHIM  is  of  saddle  stability  type,  having  (2n  —  2)-dimensional  stable  and  unstable  manifolds  IT®(/i) 
and  W^{h)  that  are  diffeomorphic  to  x  R.  Being  of  co-dimension  ^  one  with  respect  to  the  energy 

surface,  these  invariant  manifolds  act  as  separatrices,  partitioning  the  energy  surface  into  “reacting”  and 
“non-reacting”  parts  as  will  explain  in  detail  in  Sec.  2.2. 

2.1  The  normal  form 

As  mentioned  in  the  previous  section,  reaction  type  dynamics  are  induced  by  equilibrium  points  of 
saddle X Centrex  ...x centre  stability  type.  These  are  equilibria  for  which  the  matrix  associated  with  the 
linearisation  of  Hamilton’s  equations  have  eigenvalues  which  consist  of  a  pair  of  real  eigenvalues  of  equal 
magnitude  and  opposite  sign,  (-I-A,  — A),  A  S  R,  and  (n  —  I)  pairs  of  complex  conjugate  purely  imaginary 
eigenvalues,  {+iu!k,  —icok),  Wfc  €  R,  for  fc  =  2, . . . ,  n. 

The  phase  space  structures  near  equilibria  of  this  type  exist  independently  of  a  specific  coordinate  sys¬ 
tem.  However,  in  order  to  carry  out  specific  calculations  we  will  need  to  be  able  to  express  these  phase  space 
structures  in  coordinates.  This  is  where  Poincare-Birkhoff  normal  form  theory  is  used. This  is  a  well-known 
theory  and  has  been  the  subject  of  many  review  papers  and  books,  see,  e.g.,  [Deprit(1969),  Meyer(1974), 
Dragt  &  Finn(1976),  Arnol’d  et  al.(1988),  Meyer(1991),  Meyer  &  Hall(1992),  Murdock(2003)].  For  our  pur¬ 
poses  it  provides  an  algorithm  whereby  the  phase  space  structures  described  in  the  previous  section  can 

®  Briefly,  the  co-dimension  of  a  submanifold  is  the  dimension  of  the  space  in  which  the  submanifold  exists,  minus  the  dimension 
of  the  submanifold.  The  significance  of  a  submanifold  being  “co-dimension  one”  is  that  it  is  one  less  dimension  than  the  space 
in  which  it  exists.  Therefore  it  can  “divide”  the  space  and  act  as  a  separatrix,  or  barrier,  to  transport. 
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be  realised  for  a  particular  system  by  means  of  the  normal  form  transformation  which  involves  making  a 
nonlinear  symplectic  change  of  variables, 


{q,p)=TiQ,P),  (1) 

into  normal  form  coordinates,  {q,p)  which,  in  a  local  neighbourhood  C  of  the  equilibrium  point,  “decouples” 
the  dynamics  into  a  “reaction  coordinate”  and  “bath  modes.”  The  coordinate  transformation  T  is  obtained 
from  imposing  conditions  on  the  form  of  H  expressed  the  new  coordinates,  {q,p), 

Hj,F{q,p)  =  H{T-\q,p))=HiQ,P).  (2) 

These  conditions  are  chosen  such  that  Hjsif  and  the  resulting  equations  of  motions  assume  a  simple  form 
in  which  the  reaction  coordinate  and  bath  modes  “decouple” .  This  decoupling  is  one  way  of  understanding 
how  we  are  able  to  construct  the  phase  space  structures,  in  the  normal  form  coordinates,  that  govern  the 
dynamics  of  reaction. 

In  fact,  we  will  assume  that  a  (generic)  non-resonance  condition  holds  between  the  eigenvalues,  namely 
that 


k2UJ2  +  •  •  •  +  knUJn  7^  0  (3) 

for  all  integer  vectors  (fc2,...,A:„)  G  Z"“^.  When  such  a  condition  holds,  the  NF  procedure  yields  an 
explicit  expression  for  the  normalised  Hamiltonian  FInf  as  a  function  of  n  integrals  of  motion: 

H^ipiq,  p)  =  Ft'NF(F,  J2,  •  ■  •  ,  Jn)  =  \I  +  UJ2J2  UJnJn  +  h.O.t..  (4) 

Here  the  higher  order  terms  (h.o.t.)  are  at  least  of  order  two  in  the  integrals  (/,  J2, . . . ,  J„). 

The  integral,  I,  corresponds  to  a  “reaction  coordinate”  (saddle- type  DoF): 

I  =  qiPi-  (5) 

We  note  that  there  is  an  equivalent  form  of  the  reaction  coordinate:-  making  the  linear  symplectic  change  of 
variables  qi  =  {pi  —  q)/ \/2  and  pi  =  (pi  -\-  q)/ \/2,  transforms  the  above  into  the  following  form,  which  may 
be  more  familiar  to  many  readers. 


1  =  qiPi  =  ^{pI-  ql)  ■  (6) 

Geometrically  speaking,  one  can  move  freely  between  these  two  representations  by  considering  the  plane 
{qi,Pi)  and  rotating  it  by  angle  7r/4,  to  give  {q,p). 

The  integrals  for  A:  =  2, . . . ,  n,  correspond  to  “bath  modes”  (centre-type  DoF)®: 

Jk  =  ]^{pl  +  ql),  k  =  2,..., n.  (7) 

In  the  new  coordinates,  Hamilton’s  equations  have  a  particularly  simple  form: 

“^We  note  that  the  inclusion  of  diA  in  a  non-resonance  condition  would  be  vacuous;  one  cannot  have  a  resonance  of  this  kind 
between  a  real  eigenvalue,  ±A,  and  purely  imaginary  eigenvalues,  ±20;^,  fc  =  2, . . . ,  n. 

^Throughout  our  work  we  use,  somewhat  interchangeably,  terminology  from  both  chemical  reaction  dynamics  and  dynamical 
systems  theory.  This  is  most  noticable  in  our  reference  to  the  integrals  of  motion.  I  is  the  integral  related  to  reaction,  and  in 
the  context  of  dynamical  systems  theory  it  is  related  to  hyperbolic  behaviour.  The  term  “reactive  mode”  might  also  be  used 
to  describe  the  dynamics  associated  with  this  integral.  The  integrals  J25  •  •  • ,  describe  the  dynamics  associated  with  “bath 
modes” .  In  the  context  of  dynamical  systems  theory,  the  dynamics  associated  with  these  integrals  is  referred  to  as  “center  type 
dynamics”  or  “center  modes” .  A  key  point  here  is  that  integrals  of  the  motion  provide  us  with  the  natural  way  of  defining  and 
describing  the  physical  notion  of  a  “mode”.  The  nature  of  the  mode  is  defined  in  the  context  of  the  specific  application  (i.e. 
chemical  reactions)  or,  in  the  context  of  dynamical  systems  theory,  through  its  stability  properties  (i.e.  hyperbolic  or  centre). 
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9i  = 


Ai(/,  J2, . . . ,  Jn)qi, 


Pi 

ilk 

Pk 

for  k  =  2, ...  ,n,  where  we  denote 


dpi 

dqi 

dpk 

dqk 


A(/,  t/2,  .  ■  .  ,  Jn)Pl^ 

t/25  ■  ■  •  j  Jn)Pk 

^k{ky  t/2,  ■  .  ■  ,  Jn)qk} 


(8) 


A(/,  (72,  •  ■  ■  ,  (/n) 
^k{k :  J2-}  •  •  •  1  'At) 


^ArNF(/^,  J2,  ■  ■  ■  ,  Jn) 
dl 

dK^-p{I,  J2,  ■  ■  ■  ,  Jn) 

dJk 


k  =  2, . . .  ,n. 


(9) 

(10) 


The  integrals  provide  a  natural  dehnition  of  the  term  “mode”  that  is  appropriate  in  the  context  of  reaction, 
and  they  are  a  consequence  of  the  (local)  integrability  in  a  neighborhood  of  the  equilibrium  point  of  saddle- 
centre-.  .  .-centre  stability  type.  Moreover,  the  expression  of  the  normal  form  Hamiltonian  in  terms  of  the 
integrals  provides  us  a  way  to  partition  the  ‘energy’  between  the  different  modes.  We  will  provide  examples 
of  how  this  can  be  done  in  the  following. 

The  normal  form  transformation  T  in  (1)  can  be  computed  in  an  algorithmic  fashion.  One  can  give 
explicit  expression  for  the  phase  space  structures  discussed  in  the  previous  section  in  terms  of  the  normal 
form  coordinates,  {q,p).  This  way  the  phase  space  structures  can  be  constructed  in  terms  on  the  normal  form 
coordinates,  {q,p),  and  for  physical  interpretation,  transformed  back  to  the  orriginal  “physical”  coordinates, 
{Q,  P),  by  the  inverse  of  the  transformation  T.®. 

In  summary,  the  “output”  of  the  normal  form  algorithm  is  the  following: 

•  A  symplectic  transformation  T{Q,  P)  =  {q,p),  and  its  inverse  T~^(q,p)  =  (Q,  P),  that  relate  the  normal 
form  coordinates  {q,p)  to  the  original  “physical”  coordinates  {Q,P). 

•  An  expression  for  the  normalized  Hamiltonian:  in  the  form,  //nf,  in  terms  of  the  normal  form  coordi¬ 
nates  {q,p),  and  in  the  form  ATnf,  in  terms  of  the  integrals  {I,J2,...,  Jn)- 

•  Explicit  expressions  for  the  integrals  of  motion  /  and  k  =  2, . . .  ,n,  in  terms  of  the  normal  form 
coordinates. 


2.2  Explicit  definition  of  the  phase  space  structures  in  the  normal  form  coordi¬ 
nates 

As  indicated  in  the  previous  section  it  is  straightforward  to  construct  the  local  phase  space  objects  governing 
“reaction”  in  the  NF  coordinates,  {q,p).  In  this  section,  we  will  dehne  the  various  structures  in  the  normal 
form  coordinates  and  discuss  briefly  the  consequences  for  the  original  dynamical  system. 

The  structure  of  an  energy  surface  near  a  saddle  point:  For  h  <  0,  the  energy  surface  consists  of 
two  disjoint  components.  The  two  components  correspond  to  “reactants”  and  “products.”  The  top  panel 
of  Fig.  1  shows  how  the  two  components  project  to  the  various  planes  of  the  normal  form  coordinates.  The 
projection  to  the  plane  of  the  saddle  coordinates  {qi,pi)  is  bounded  away  from  the  origin  by  the  two  branches 
of  the  hyperbola,  qipi  =  /  <  0,  where  /  is  given  implicitly  by  the  energy  equation  with  the  centre  actions 
Jfc,  k  =  2, . . . ,  n,  set  equal  to  zero:  ATnf)/,  0, . . . ,  0)  =  h  <  0.  The  projections  to  the  planes  of  the  centre 
coordinates,  {qk,Pk),  k  =  2, . . .  ,n,  are  unbounded. 

®The  original  coordinates  (Q,  P)  typically  have  an  interpretation  as  configuration  space  coordinates  and  momentum  coordi¬ 
nates.  The  normal  form  coordinates  {q,p),  in  general,  do  not  have  such  a  physical  interpretation  since  both  q  and  p  are  nonlinear 
functions  of  both  Q  and  P. 


6 


P2A 


^2 


P2A 


P3A 


^2 


P2  A 


P3  A 


^12 


% 


^3 


Figure  ll  Projection  of  energy  surfaces  (turquoise  regions)  to  the  planes  of  the  normal  form  coordinates.  The  energy  surface 
have  energy  h  <  0  (top  panel),  h  =  0  (middle  panel),  h  >  0  (bottom  panel)  . 
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Figure  2:  Projection  of  the  dividing  surface  and  reacting  and  nonreacting  trajectories  to  the  planes  of  the  normal  form 
coordinates.  In  the  plane  of  the  saddle  coordinates,  the  projection  of  the  dividing  surface  is  the  dark  red  diagonal  line  segment, 
which  has  qi  =  pi-  In  the  planes  of  the  centre  coordinates,  the  projections  of  the  dividing  surface  are  the  dark  red  discs. 
Forward  and  backward  reactive  trajectories  (yellow  and  blue)  project  to  the  first  and  third  quadrant  in  the  plane  of  the  saddle 
coordinates,  respectively,  and  pass  through  the  dividing  surface.  The  red  and  green  curves  mark  nonreactive  trajectories  on  the 
reactant  side  (pi  —  qi  >0),  and  on  the  product  side  (pi  —  qi  <  0),  of  the  dividing  surface,  respectively.  The  turquoise  regions 
indicate  the  projections  of  the  energy  surface. 


At  =  0,  the  formerly  disconnected  components  merge  (the  energy  surface  bifurcates),  and  for  h  >  0^ 
the  energy  surface  has  locally  the  structure  of  a  spherical  cylinder,  x  R.  Its  projection  to  the  plane 

of  the  saddle  coordinates  now  includes  the  origin.  In  the  first  and  third  quadrants  it  is  bounded  by  the  two 
branches  of  the  hyperbola,  qipi  =  /  >  0,  where  I  is  again  given  implicitly  by  the  energy  equation  with 
all  centre  actions  equal  to  zero,  but  now  with  an  energy  greater  than  0:  0, . . . ,  0)  =  h  >  0.  The 

projections  to  the  planes  of  the  centre  coordinates  are  again  unbounded.  This  is  illustrated  in  the  bottom 
panel  of  Fig.  1. 

The  dividing  surface,  and  reacting  and  nonreacting  trajectories:  On  an  energy  surface  with  h  >  0, 
we  define  the  dividing  surface  by  gi  =  pi-  This  gives  a  (2n  —  2)-sphere  which  we  denote  by  Its 

projection  to  the  saddle  coordinates  simply  gives  a  line  segment  through  the  origin  which  joins  the  boundaries 
of  the  projection  of  the  energy  surface,  as  shown  in  Fig.  2.  The  projections  of  the  dividing  surface  to  the 
planes  of  the  centre  coordinates  are  bounded  by  circles  (pl+Qk)/^  —  k  =  2, . . . ,  d,  where  Jk  is  determined 
by  the  energy  equation  with  the  other  centre  actions,  Ji,  I  ^  k,  and  the  saddle  integral,  7,  set  equal  to  zero. 
The  dividing  surface  divides  the  energy  surface  into  two  halves,  pi  —  <71  >  0  and  pi  —  qi  <  0,  corresponding 
to  reactants  and  products. 

As  mentioned  above,  trajectories  project  to  hyperbolae  in  the  plane  of  the  saddle  coordinates,  and  to 
circles  in  the  planes  of  the  centre  coordinates.  The  sign  of  I  determines  whether  a  trajectory  is  nonreacting 
or  reacting,  see  Fig.  2.  Trajectories  which  have  /  <  0  are  nonreactive  and  for  one  branch  of  the  hyperbola 
qiPi  =  I  they  stay  on  the  reactants  side  and  for  the  other  branch  they  stay  on  the  products  side;  trajectories 
with  /  >  0  are  reactive,  and  for  one  branch  of  the  hyperbola  qiPi  =  I  they  react  in  the  forward  direction, 
i.e.,  from  reactants  to  products,  and  for  the  other  branch  they  react  in  the  backward  direction,  i.e.,  from 
products  to  reactants.  The  projections  of  reactive  trajectories  to  the  planes  of  the  centre  coordinates  are 
always  contained  in  the  projections  of  the  dividing  surface.  In  this,  and  other  ways,  the  geometry  of  the 
reaction  is  highly  constrained.  There  is  no  analogous  restriction  on  the  projections  of  nonreactive  trajectories 
to  the  centre  coordinates. 

The  normally  hyperbolic  invariant  manifold  (NHIM)  and  its  relation  to  the  ‘activated  complex’: 

On  an  energy  surface  with  h  >  0,  the  NHIM  is  given  by  qi  =  pi  —  0.  The  NHIM  has  the  structure  of  a 
(2n  —  3)-sphere,  which  we  denote  by  §NfflM(^)-  ^he  NHIM  is  the  equator  of  the  dividing  surface;  it  divides 
it  into  two  “hemispheres” :  the  forward  dividing  surface,  which  has  qi  =  pi  >  0,  and  the  backward  dividing 
surface,  which  has  qi  =  pi  <  0.  The  forward  and  backward  dividing  surfaces  have  the  structure  of  (2n  —  2)- 
dimensional  balls,  which  we  denote  by  and  B‘^~^(h),  respectively.  All  forward  reactive  trajectories 

cross  backward  reactive  trajectories  cross  B'^f^ih).  Since  qi  =  pi  —  0  in  the  equations  of 

motion  (8)  implies  that  qi  =  pi  =  0,  the  NHIM  is  an  invariant  manifold,  i.e.,  trajectories  started  in  the 
NHIM  stay  in  the  NHIM  for  all  time.  The  system  resulting  from  gi  =  pi  =  0  is  an  invariant  subsystem  with 
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Figure  3:  The  projection  of  the  NHIM  and  the  local  parts  of  its  stable  and  unstable  manifolds,  and  to  the 

planes  of  the  normal  form  coordinates.  In  the  plane  of  the  saddle  coordinates,  the  projection  of  the  NHIM  is  the  origin  marked 
by  the  blue  bold  point,  and  the  projection  of  and  are  the  pi-axis  and  gi-axis,  respectively.  W^{h)  consists  of  the 

forward  and  backward  branches  Wj  (/i)  and  W^{h),  which  have  pi  >  0  and  pi  <  0,  respectively;  W'^{h)  consists  of  W^{h)  and 
W^{h),  which  have  gi  >  0  and  qi  <  0,  respectively.  In  the  plane  of  the  centre  coordinates,  the  projections  of  the  NHIM,  iy®(/i), 
and  iy“(/i)  (the  blue  circular  discs)  coincide  with  the  projection  of  the  dividing  surface  in  Fig.  2.  The  turquoise  regions  mark 
the  projections  of  the  energy  surface. 


one  degree  of  freedom  less  than  the  full  system.  In  fact,  qi  =  pi  =  0  defines  the  centre  manifold  associated 
with  the  saddle-centre-’  •  •  -centre  equilibrium  point,  and  the  NHIM  at  an  energy  h  greater  than  the  energy  of 
the  quilibrium  point  is  given  by  the  intersection  of  the  centre  manifold  with  the  energy  surface  of  this  energy 
h  [Uzer  et  al.(2002),  Waalkens  &  Wiggins(2004)]. 

This  invariant  subsystem  with  one  degree  of  freedom  less  than  the  full  system  is  the  “activated  complex” 
(in  phase  space),  located  between  reactants  and  products.  The  NHIM  can  be  considered  to  be  the  energy 
surface  of  the  activated  complex.  In  particular,  all  trajectories  in  the  NHIM  have  1  =  0. 

The  equations  of  motion  (8)  also  show  that  pi  —  ()i  <  0  on  the  forward  dividing  surface  and 

Pi  —  qi  >  0  on  the  backward  dividing  surface  Hence,  except  for  the  NHIM,  which  is  is  an  invariant 

manifold,  the  dividing  surface  is  everywhere  transverse  to  the  Hamiltonian  flow.  This  means  that  a  trajectory, 
after  having  crossed  the  forward  or  backward  dividing  surface,  or  B'^~^{h),  respectively,  must  leave 

the  neighbourhood  of  the  dividing  surface  before  it  can  possibly  cross  it  again.  Indeed,  such  a  trajectory 
must  leave  the  local  region  in  which  the  normal  form  is  valid  before  it  can  possibly  cross  the  dividing  surface 
again. 

The  NHIM  has  a  special  structure:  due  to  the  conservation  of  the  centre  actions,  it  is  filled,  or  foliated, 
by  invariant  (n  —  l)-dimensional  tori,  T”“^.  More  precisely,  for  d  =  3  degrees  of  freedom,  each  value  of  J2 
implicitly  defines  a  value  of  J3  by  the  energy  equation  itrNF(0,  J2,Jz)  =  h.  For  three  degrees  of  freedom,  the 
NHIM  is  thus  foliated  by  a  one-parameter  family  of  invariant  2-tori.  The  end  points  of  the  parameterization 
interval  correspond  to  J2  =  0  (implying  q2  =  P2  =  0)  and  J3  =  0  (implying  qs  =  ps  =  0),  respectively.  At 
the  end  points,  the  2-tori  thus  degenerate  to  periodic  orbits,  the  so-called  Lyapunov  periodic  orbits. 

The  stable  and  unstable  manifolds  of  the  NHIM  forming  the  phase  space  conduits  for  reactions: 

Since  the  NHIM  is  of  saddle  stability  type,  it  has  stable  and  unstable  manifolds,  W^{h)  and  W'^{h).  The 
stable  and  unstable  manifolds  have  the  structure  of  spherical  cylinders,  x  R.  Each  of  them  consists  of 

two  branches:  the  “forward  branches”,  which  we  denote  by  Wj{h)  and  IF“(/i),  and  the  “backward  branches”, 
which  we  denote  by  VF^(h)  and  Wff{h).  In  terms  of  the  normal  form  coordinates,  W({h)  is  given  by  <71  =  0 
with  Pi  >  0,  Wf{h)  is  given  by  pi  =0  with  qi  >  0,  W§{h)  is  given  by  <71  =  0  with  pi  <  0,  and  W{f{h)  is 
given  by  Pi  =  0  with  qi  <  0,  see  Fig.  3.  Trajectories  on  these  manifolds  have  1  =  0. 

Since  the  stable  and  unstable  manifolds  of  the  NHIM  are  of  one  less  dimension  than  the  energy  surface, 
they  enclose  volumes  of  the  energy  surface.  We  call  the  union  of  the  forward  branches,  Wj{h)  and  Wf“(h), 
the  forward  reactive  spherical  cylinder  and  denote  it  by  Wi(h).  Similarly,  we  define  the  backward  reactive 
spherical  cylinder,  Wh{h),  as  the  union  of  the  backward  branches,  Wf^{h)  and  W^{h). 

The  reactive  volumes  enclosed  by  Wf{h)  and  Wb{h)  are  shown  in  Fig.  4  as  their  projections  to  the  normal 
form  coordinate  planes.  In  the  plane  of  the  saddle  coordinates,  the  reactive  volume  enclosed  by  Wf{h)  projects 
to  the  first  quadrant.  This  projection  is  bounded  by  the  corresponding  hyperbola  qipi  =  I,  with  I  obtained 
from  Kj<iF{I,0, . . .  ,0)  =  h.  Likewise,  Wb{h)  projects  to  the  third  quadrant  in  the  ((7i,pi)-plane.  Wf{h) 
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Figure  4:  Projections  of  the  reactive  volumes  enclosed  by  the  forward  and  backward  reactive  spherical  cylinders,  Wf{h)  and 
s,nd  the  forward  and  backward  reactions  paths,  to  the  planes  of  the  normal  form  coordinates.  The  volumes  enclosed  by 
Wf(h)  and  Wb(h)  project  to  the  dark  pink  and  green  regions  in  the  first  and  third  quadrant  in  the  plane  of  the  saddle  coordinates, 
respectively.  These  volumes  project  to  the  dark  green/dark  pink  brindled  disks  in  the  planes  of  the  centre  coordinates,  where 
their  projections  coincide  with  the  projection  of  the  NHIM  and  the  dividing  surface  in  Figs.  2  and  3.  The  forward  and  backward 
reaction  paths  project  to  the  two  branches  of  a  hyperbola  marked  blue  in  the  first  and  third  quadrant  in  the  plane  of  the  saddle 
coordinates,  respectively,  and  to  the  origins  (bold  blue  points)  in  the  planes  of  the  centre  coordinates.  The  turquoise  regions 
mark  the  projections  of  the  energy  surface. 


encloses  all  forward  reactive  trajectories;  Wb{h)  encloses  all  backward  reactive  trajectories.  All  nonreactive 
trajectories  are  contained  in  the  complement. 

Forward  and  backward  reaction  paths:  The  local  geometry  of  Wf{h)  and  Wh{h)  suggests  a  natural 
definition  of  dynamical  forward  and  backward  reaction  paths  as  the  unique  paths  in  phase  space  obtained  by 
putting  all  of  the  energy  of  a  reacting  trajectory  into  the  reacting  mode,  i.e.,  setting  <72  =  •  •  •  =  =  P2  =  •  •  •  = 

Pn  =  0.  This  gives  the  two  branches  of  the  hyperbola  qipi  =  I,  with  I  obtained  from  0, . . . ,  0)  =  /i, 

which  in  phase  space  are  contained  in  the  plane  of  the  saddle  coordinates,  see  Fig.  4.  This  way,  the  forward 
(respectively,  backward)  reaction  path  can  be  thought  of  as  the  “centre  curve”  of  the  relevant  volume  enclosed 
by  the  forward  (resp.,  backward)  reactive  spherical  cylinder  Wf{h)  (resp.,  Wh{h)).  These  reaction  paths  are 
the  special  reactive  trajectories  which  intersect  the  dividing  surface  at  the  “poles”  (in  the  sense  of  North  and 
South  poles,  where  qi  =  pi  assumes  its  maximum  and  minimum  value  on  the  dividing  surface). 

The  Transmission  Time  through  the  Transition  State  Region:  The  normal  form  coordinates  provide 
a  way  of  computing  the  time  for  all  trajectories  to  cross  the  transition  region.  We  illustrate  this  with  a  forward 
reacting  trajectory  (a  similar  argument  and  calculation  can  be  applied  to  backward  reacting  trajectories). 
We  choose  the  boundary  for  the  entrance  to  the  reaction  region  to  be  pi  —  gi  =  c  for  some  constant  c  >  0, 
i.e.,  initial  conditions  which  lie  on  the  reactant  side  of  the  transition  state,  and  the  boundary  for  exiting  the 
reaction  region  to  be  pi  —  qi  =  —c  on  the  product  side.  We  now  compute  the  time  of  flight  for  a  forward 
reacting  trajectory  with  initial  condition  on  pi  —  qi  =  c  to  reach  pi  —  qi  =  —c  on  the  product  side.  The 
solutions  are  qi{t)  =  qi{0)  exp(A(/,  J2,. . . ,  Jn)t)  and  pi{t)  =  pi(0)  exp(— A(/,  J2, . . . ,  Jn)t)  (see  (8)),  where 
A(/,  J2,  •  •  ■ ,  Jn)  is  determined  by  the  initial  conditions.  This  gives  the  time  of  flight  as 

r=(A(/,J2,...yn))-'ln(^^)  .  (11) 

The  time  diverges  logarithmically  as  qi(0)  ^  0,  i.e.,  the  closer  the  trajectory  starts  to  the  boundary  Wi{h). 
It  is  not  difficult  to  see  that  the  time  of  flight  is  shortest  for  the  centre  curve  of  the  volume  enclosed  by  Wi{h), 
i.e.,  the  trajectory  which  traverses  the  transition  state  region  fastest  is  precisely  our  forward  reaction  path.  A 
similar  construction  applies  to  backward  reactive  trajectories. 

In  fact,  the  normal  form  can  be  used  to  map  trajectories  through  the  transition  state  region,  i.e.  the 
phase  space  point  at  which  a  trajectory  enters  the  transition  state  region  can  be  mapped  analytically  to  the 
phase  space  point  at  which  the  trajectories  exits  the  transition  state  region. 
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2.3  The  foliation  of  the  reaction  region  by  Lagrangian  submanifolds 

The  existence  of  the  n  integrals  of  motion,  {I,  J2,  ■  ■  ■ ,  Jn),  induce  phase  space  structures  which  lead  to  further 
constraints  on  the  trajectories  in  addition  to  the  ones  described  above.  In  order  to  describe  these  structures 
and  the  resulting  constraints  it  is  useful  to  introduce  the  so  called  momentum  map,  M  [Guillemin(1994), 
Marsden  &  Ratiu(1999)]  which  maps  a  point  (gi, . . . ,  qn,Pi,  ■  ■  ■  ,Pn)  in  the  phase  space  R”  x  R”  to  the  n 
integrals  evaluated  at  this  point: 

M:R"xR"^R",  (gi, . . . ,  . . .  ,p„)  ^  (7,  J2,  •  ■  • ,  •  (12) 

The  preimage  of  a  value  for  the  constants  of  motion  {I,  J2, .  ■ . ,  Jn)  under  M  is  called  a  fibre.  A  fibre  thus 
corresponds  to  the  common  level  set  of  the  integrals  in  phase  space. 

A  point  {qi, . . .  ,qn,pi,  ■ .  ■  ,Pn)  is  called  a  regular  point  of  the  momentum  map  if  the  linearisation  of  the 
momentum  map,  DA4,  has  full  rank  n  at  this  point,  i.e.,  if  the  gradients  of  the  n  integrals  I,  Jfe,  k  =  2, . . .  ,n, 
with  respect  to  the  phase  space  coordinates  {q,p)  are  linearly  independent  at  this  point.  If  the  rank  of  DAI 
is  less  than  n  then  the  point  is  called  irregular.  A  hbre  is  called  regular  if  it  consists  of  regular  points  only. 
Else  it  is  called  an  irregular  fibre.  In  fact  almost  all  fibres  are  regular.  They  are  n-dimensional  manifolds 
given  by  the  Cartesian  product  of  an  hyperbola  qipi  =  /  in  the  saddle  plane  {qi,pi)  and  n  —  1  circles  in 
the  centre  planes  {qk,Pk)i  k  =  2, . . .  ,n.  Since  the  hyperbola  qiPi  =  I  consists  of  two  branches  each  of  which 
have  the  topology  of  a  line  R,  the  regular  fibres  consist  of  two  disjoint  toroidal  cylinders,  x  R,  which 

are  the  Cartesian  products  of  a  (n  —  l)-dimensional  torus  and  a  line.  We  denote  these  toroidal  cylinders  by 

^  \{P2  +  (ll)  =  ,^{pl  +  ql)  =  Jn,qi  >  0}  (13) 

and 

=  {i9,p)  G  18.^”  :  pm  =  i,^{pl  +  ql)  =  'h,---  ,^{pl  +  ql)  =  J„,gi  <  o}.  (14) 

and  Aj  are  Lagrangian  manifolds  [Arnold(1978)].  The  Lagrangian  manifolds  consists  of  all 

trajectories  which  have  the  same  constants  of  motion.  In  particular  the  Lagrangian  manifolds  are  invariant, 
i.e.  a  trajectory  with  initial  condition  on  a  Lagrangian  manifold  will  stay  in  the  Lagrangian  manifold  for  all 
time.  For  /  <  0,  the  Lagrangian  manifolds  j  and  A^j^  j  consist  of  nonreactive  trajectories  in  the 

reactants  resp.  products  components  of  the  energy  surface.  For  /  >  0,  A^j^  consists  of  forward  reactive 
trajectories,  and  Aj"j^  j  consists  of  backward  reactive  trajectories. 

We  will  discuss  the  irregular  hbres  in  more  detail  in  Sec.  4. 

2.4  Implications  for  the  original  system 

The  normalization  procedure  proceeds  via  formal  power  series  manipulations  whose  input  is  a  Taylor  expan¬ 
sion  of  the  original  Hamiltonian,  H ,  necessarily  up  to  some  finite  order,  M,  in  homogeneous  polynomials. 
For  a  particular  application,  this  procedure  naturally  necessitates  a  suitable  choice  of  the  order,  M ,  for  the 
normalization,  after  which  one  must  make  a  restriction  to  some  local  region,  C,  about  the  equilibrium  point  in 
which  the  resulting  computations  achieve  some  desired  accuracy.  Hence,  the  accuracy  of  the  normal  form  as  a 
power  series  expansion  truncated  at  order  M  in  a  neighborhood  C  is  determined  by  comparing  the  dynamics 
associated  with  the  normal  form  with  the  dynamics  of  the  original  system.  There  are  several  independent 
tests  that  can  be  carried  out  to  verify  accuracy  of  the  normal  form.  Straightforward  tests  that  we  use  are 
the  following: 

•  Examine  how  the  integrals  associated  with  the  normal  form  change  on  trajectories  of  the  full  Hamilto¬ 
nian  (the  integrals  will  be  constant  on  trajectories  of  the  normal  form). 

•  Check  invariance  of  the  different  invariant  manifolds  (i.e.  the  NHIM  and  its  stable  and  unstable 
manifolds)  with  respect  to  trajectories  of  the  full  Hamiltonian. 

Both  of  these  tests  will  require  us  to  use  the  transformations  between  the  original  coordinates  and  the 
normal  form  coordinates.  Specihc  examples  where  M,  C  and  accuracy  of  the  normal  forms  are  considered  can 
be  found  in  [Waalkens  et  al. (2004a),  Waalkens  et  al. (2004b),  Waalkens  et  al. (2005a),  Waalkens  et  al. (2005b), 
Waalkens  et  al.  (2005c)]. 
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3  Models  of  the  local  phase  space  structures  for  2-DoF  systems 


In  the  discussion  of  the  phase  space  structures  governing  reaction  type  dynamics  in  Sec.  2.2  we  visualized 
these  structures  by  projecting  them  to  the  2-dimensional  planes  of  the  coordinates  and  conjugate  momenta 
of  the  normal  form.  The  topology  of  the  phase  space  structures  and  the  interrelationship  between  them  is 
to  a  large  extent  obscured  in  these  projection  because  of  the  lacking  injectivity  of  the  projections.  For  n  =  2 
DoFhowever,  a  surface  of  constant  energy  containing  the  phase  space  structures  is  three-dimensional.  For 
this  case,  it  is  therefore  possible  to  visualize  the  phase  space  structures  in  three-dimensional  space.  In  this 
section  we  will  develop  three  models  to  explicitly  implement  such  visualizations.  Our  discussion  will  be  in 
the  context  of  linear  Hamiltonian  systems  since  they  allow  us  to  describe  the  phase  space  structures  in  terms 
of  explicit  formulae.  The  influence  of  nonlinear  terms  is  discussed  in  Sec.  3.4. 

Therefore  we  begin  by  considering  the  quadratic  Hamiltonian 


H{x,y,p^,Py)  = 

with  A,  w  >  0. 

This  gives  the  linear  Hamiltonian  vector  field 

X  = 

Px  = 

y  = 

Py  = 


(15) 

dH  , 

^Px 

(16) 

dH  , 

=  Ax, 
ox 

(17) 

dH 

(18) 

II 

1 

e 

(19) 

The  Hamiltonian  (15)  is  already  in  normal  form.  To  keep  the  notation  simple  we  will  not  specifically 
indicate  this  by  a  subscript  NF.  In  fact  we  have 

H{x,y,Px,Py)  =  K{I,J)  =  XI +  ujJ ,  (20) 

where  /  and  J  are  the  integrals  given  in  (6))  and  (7),  respectively. 

For  a  fixed  energy  h,  the  three  dimensional  energy  surface  is  given  by  the  level  set  of  the  Hamiltonian  (15) 
as  follows: 


j:{h)  =  {ix,y,p^,Py)  G  :  \{pI  -  +uj{pl+y'^)  =  2h}  .  (21) 

We  rewrite  the  energy  equation  in  the  form 

Xpl  +  Lopy  +  ojy"^  =  2h  +  Xx^  ,  (22) 

and  fix  x  such  that  the  right  hand  side  of  this  equation  is  positive.  Note  that  if  h  >  0  the  right  hand  side 
will  be  positive  for  any  fixed  x  G  R,  while  if  <  0  it  will  be  positive  only  for  |x|  >  yj —2h/X,  i.e.  for  x 
sufficiently  small  or  large.  For  such  an  x,  Eq.  (22)  defines  an  ellipsoid  in  the  three-dimensional  Euclidean 
space,  R^,  with  coordinates  {px,Py,y)-  This  ellipsoid  is  rotationally  symmetric  about  the  Pa;-axis.  It  has  two 

semi-axes  of  length  ((2/i  -|-  Ax^)/a;)^^^  and  one  semi-axis  of  length  ((2/i  -|-  Ax^)/A)^^^.  Topologically,  this  is 
a  two-dimensional  sphere,  which  we  denote  by  §^(/i),  i.e. 

=  {{Px,Py,y)  &R^  ■■  Xpl+u}{pl  +  y‘^)  =2h  + Xx"^}  .  (23) 

The  energy  surface  thus  appears  to  have  the  structure  of  a  one-parameter  family  of  two-dimensional  spheres 
parametrized  by  x,  where  the  ‘radius’  (the  length  of  the  shortest  semi-axis  of  the  ellipsoid  (22))  can  shrink  to 
zero  if  /i  <  0.  The  three  models  of  the  energy  surface  and  the  phase  space  structures  contained  in  it  that  we 
present  in  the  following  differ  by  the  way  of  representing  the  2-spheres  For  clarity,  let  us  summarize 
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the  phase  space  structures  discussed  in  Sec.  2  for  the  Hamiltonian  (15)  before  we  proceed.  We  will  restrict 
ourselves  to  the  case  h  >  0. 

The  dividing  surface  is  obtained  by  setting  a;  =  0  on  the  energy  surface.  This  gives  the  two-dimensional 
sphere 


^dsih)  =  S(h)  n  {a;  =  0}  =  {{x,y,Pa;,Py)  G  :  \pl  +  uj {pi  +  y'^)  =  2h,  x  =  0}  .  (24) 

For  the  2-DoF  case,  the  NHIM  is  an  unstable  periodic  orbit  (the  Lyapunov  periodic  orbit  associated  with  the 
saddle  point).  For  the  Hamiltonian  (15),  it  is  obtained  by  setting  a;  =  0  and  Px  =  0  on  the  energy  surface. 
This  gives  the  circle 


^nhimW  =  n  {x  =  0,  Pa:  =  0}  =  {{x,y,Px,Py)  G  :  uj{pl+y‘^)  =  2h,  x  =  0,  =  0}  .  (25) 

The  NHIM  divides  the  dividing  surface,  S^g(/i),  into  the  forward  hemisphere 

=  S(h)  n  {x  =  0,  p,^  >  0}  =  {ix,y,Px,Py)  G  R"*  :  Px  =  \J {2h  -  dj{pl  +  y"^)) / X,  x  =  0}  (26) 

and  the  backward  hemisphere 

=  S(h)  n  {x  =  0,  Pj,  <  0}  =  {{x,y,Px,Py)  G  R"*  :  Px  =  - {2h  -  uj{pl  +  y"^)) / X,  x  =  0}  ,  (27) 

which  are  two  dimensional  balls  or  disks.  The  stable  manifold  IT®  of  Snhim(^)  consists  of  the  forward  branch 

Wf®  =  S(h)  n  {x  =  -px,  X  <  0}  =  {(x,p,px,py)  G  R^  :  uj{pl  +  y^)  =  2h,  x  =  -px,  x  <  0}  (28) 

and  the  backward  branch 

=  T,{h)  n  {x  =  -px,  X  >  0}  =  {{x,y,Px,Py)  G  R^  :  uj{pl  +  =  2h,  x  =  -px,  x  >  0}  .  (29) 

Similarly,  the  unstable  manifold  1T“  of  SJjhim  (^)  consists  of  the  forward  branch 

Wf"  =  S(h)  n  {x  =pa;,  X  >  0}  =  {{x,y,Px,Py)  G  R'^  :  uj{pl  +  y'^)  =  2h,  x  =  Px,  x  >  0}  (30) 

and  the  backward  branch 

VFb  =  S(h)  n  {x  =  Pa;,  X  <  0}  =  {{x,y,Px,Py)  G  R^  :  uj{pl  +  y^)  =  2h,  x  =  Px,  x  <  0}  .  (31) 

Finally,  the  forward  and  backward  reaction  paths  are  obtained  from  putting  all  the  energy  in  the  x  DoF. 
This  gives  the  lines 

S(/i)  n  {y  =  py  =  0,  Px  >  0}  =  {{x,y,Px,Py)  G  R'^  :  Px  =  \/(2h -f  Ax)/A,  y  =  pj^  =  0}  (32) 

for  the  forward  reaction  path,  and 

Y,{h)  n  {y  =  Py  =  0,  Px  <  0}  =  {(x,y,Px,Py)  G  R'^  :  Px  =  -a/ {2h  +  Ax)/A,  y  =  pj^  =  0}  (33) 

for  the  backward  reaction  path. 

Note  that  the  energy  surface  (21),  the  phase  space  structures  (24)-(33)  and  also  the  Lagrangian  cylinders 
defined  in  (13)  and  (14)  are  invariant  under  the  symmetry  action 

(x,y,px,Pj/)  ^  (x,e"*‘^y,Px,e"*>j;) ,  p  G  (34) 
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3.1  Model  I:  projection  of  2-spheres  to  their  equatorial  planes 

The  main  idea  of  the  first  model  is  to  view  the  2-sphere  (23)  as  consisting  of  two  hemispheres  and  an  equator, 
i.e.  as  the  disjoint  union  of  two  open  two-dimensional  disks  or  balls,  B^,  and  a  circle,  or  equivalently  one¬ 
dimensional  sphere,  S^,  and  using  parametrizations  of  these  geometric  structures  to  parametrize  the  energy 
surface. 

The  decomposition  of  a  2-sphere,  S^,  in  terms  of  2-balls  and  a  circle  is  familiar  from  visualizing  the  surface 
of  the  Earth  by  drawing  two  filled-in  disks  one  each  for  the  Northern  and  Southern  hemispheres,  and 

identifying  the  boundaries  of  the  two  disks, 

=  -^North  U  ^Equator  ^  -Ssouthi  (35) 

which  we  illustrate  in  Fig.  5. 


Figure  5:  A  model  of  the  sphere  formed  by  taking  two  open  discs  (2-balls)  i?North  -^South  identifying 
points  on  their  boundaries;  the  “equator”  SEquator-  ^1®°  shown  is  the  representation  of  a  smooth  closed  curve 
on  which  passes  through  the  North  and  South  poles,  N  and  S';  the  points  at  which  the  representation  of 
the  curve  “jumps”  from  one  disc  in  the  model  to  the  other  disc  are  marked  by  dots. 


For  the  2-spheres  S^(h)  in  Eq.  (23),  we  define  the  decomposition  (35)  as 

where 


(36) 


^+,x(^)  =  {{Px,Py,y)^'R^  ■  +  ^ {2h  +  Ax^  -  ujpl  -  ujy"^) / \  >  0}  ,  (37) 

=  {{Px,Py,y)  G  -f  2/^)  <  2/i-  Aa;^  {2h  +  Xx"^  -  ujpl  -  ujy"^) / X  <  0}  ,  (38) 

^l,x{^)  =  {(Pa;,Py,y)  G  -f  2/2)  =  2/i- Aa;^  Pa,  =  0}  .  (39) 

We  can  view  the  manifolds  B^  ^{h)  and  B^  ^{h)  as  graphs  of  the  functions 

"Pi.x  : -B^(/i)  ^  B.,  {y,Py)  ^±^{‘2h  +  Aa;2  -  wpg  -  ujy'^)/\  ,  (40) 

where  the  domains  are  the  circular  disks 

Bl{h)  =  {{Py,y)  G  R2  :  +  <  2h  +  \x'^}  C  .  (41) 

The  graphs  ‘join’  at  the  equator  Sj  „,(/i).  In  order  to  geometrically  visualize  S2(/i)  we  take  two  copies  of  the 
circular  disks  B^(h)  in  the  {y,py)-plane  where  p^  is  given  by  the  function  V+^x  on  one  copy  which  we  denote 
by  B^  xW  s-M  the  function  V-^x  on  the  other  copy  which  we  denote  by  ^^(h).  B^^^{h)  thus  denote  the 
images  of  B\^^{h)  under  the  projection 

'^P^.x  ■■  B\^^{h)  ,  {px,Py,y)  ^  {y,Py) ,  (42) 
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Figure  6:  Model  I:  projection  of  to  the  (y,py)-plane. 

see  Fig.  6.  Along  the  boundaries  of  and  we  have  Px  =  0,  and  the  points  on  the  boundaries 

which  have  the  same  {y,py)  are  identified. 

The  first  model  (Model  I)  of  the  energy  surface  and  the  phase  space  structures  contained  in  it  is  obtained 
from  defining  a  map  from  the  energy  surface  S(/i)  to  using  the  projections  'np^x,  x  G  R,  in  (42)  to  map 
points  {x,y,Px,Py)  G  ^{h)  according  to 


:  {x,y,Px,Py)  {x,TTp^x{Px,Pv,y))  =  {x,y,py) .  (43) 

The  map  is  singular  along  the  cylinder  given  by  the  union  of  the  circles  ^.(ft.)  in  (39)  with  x  G  R.  This 
singular  cylinder  is  mapped  to  the  quadric 


2h/u}  2h/X  2h/uj 


(44) 


which  is  a  single  sheeted  hyperboloid  that  is  rotationally  symmetric  about  the  x  axis.  ^  Each  point  in 
(x,  y,Py)  G  R^  with 


2h/uj  2h/X  2h/uj 


<  1 


(45) 


has  two  preimages  that  differ  by  the  sign  of  Px-  We  can  thus  represent  the  energy  surface  E(ft,)  in  the  three 
dimensional  space,  R^,  in  terms  of  the  manifolds 


E^±(/i)  =  {{x,y,py)  G  R^ 
=  {ix,y,Py)  G  R^ 


X  G  R,  {y,py)  G  Bl^xih)} 


pI 


2h/oj  2h/X  2h/uj 


<1}. 


(46) 


where  Pa;  =  (2ft,  +  x^— is  positive  on  E5|_(/i)  andpa;  =  — (2/i  +  x^— p^  — is  negative  on  SL(/i). 
The  full  energy  surface  E(ft,)  is  obtained  by  gluing  S^(/i)  along  their  boundaries  (44)  where  Px  =  0.  We  show 
illustrations  of  E5|_(ft,)  and  Y}_{h)  in  the  top  panels  of  Fig.  7.  As  a  consequence  of  the  symmetry  (34)  the 
components  S^(/i)  are  both  rotationally  symmetric  about  the  x-axis.  It  is  therefore  useful  to  consider  also 
the  intersections  of  E5|_(/i)  and  E5|_(/i)  (and  of  the  other  phase  space  structures  discussed  in  the  following)  with 
the  plane  y  =  0  as  shown  in  the  bottom  panels  of  Fig.  7.  It  is  easy  to  see  the  wide-narrow-wide  geometry  of 
the  bottleneck  type  structure  of  the  energy  surface  in  this  representation  of  the  energy  surface.  The  reactants 
region  is  on  the  left,  the  product  regions  is  on  the  right. 

^We  note  that  for  h  <  0  the  quadrics  (44)  are  two-sheeted  hyperboloids.  For  h  =  0,  the  quadric  is  a  ‘diabolo’. 
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Figure  7:  Model  I:  the  energy  surface.  The  orientation  of  the  x,  y  and  py  axes  is  indicated  by  the  dreibein. 
(Note  that  the  origin  of  the  coordinate  system  is  not  the  center  of  the  dreibein.)  The  top  left  panel  shows 
S5i_(/i);  the  top  left  panel  shows  SL(/i).  The  two  panels  at  the  bottom  show  the  intersections  of  S5|_(/i)  and 
Y}_(h)  with  the  plane  y  =  0. 
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Figure  8:  Model  I:  nonreactive  trajectories  that  stay  in  the  reactants  component.  The  red  wireframes  in  the 
top  panels  indicate  the  two  pieces  of  a  Lagrangian  cylinder  Aj  j  defined  in  (14)  for  J  <  0  contained  in 
and  SL(/i),  respectively.  The  black  line  shows  one  of  the  trajectories  contained  in  Aj  j.  The  red  lines  in  the 
bottom  panels  show  the  intersection  of  Aj  j  with  the  plane  y  =  0. 
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Figure  9:  Model  I:  nonreactive  trajectories  that  stay  in  the  products  component.  The  green  wireframes  in 
the  top  panels  indicate  the  two  pieces  of  a  Lagrangian  cylinder  defined  in  (13)  for  /  <  0  contained  in 
and  IlL(h),  respectively.  The  black  line  shows  one  of  the  trajectories  contained  in  Af  j.  The  green 
lines  in  the  bottom  panels  show  the  intersection  of  with  the  plane  y  =  0. 


In  Fig.  8  we  show  the  image  under  of  a  Lagrangian  cylinder  Ajj  defined  in  (14)  for  /  <  0  together 
with  a  nonreactive  trajectory  contained  in  this  cylinder.  Note  that  for  a  2-DoF  system  and  given  values  for 
the  integral  /  and  the  energy  h,  the  second  integral  J  is  fixed  by  the  energy  equation  K{I,  J)  =  h  (see  (20)). 
The  trajectory  approaches  the  bottleneck  region  with  positive  px  from  the  reactants  component  of  the  energy 
surface  within  S?|_(/i)  (the  left  panels  of  Fig.  8),  reaches  the  boundary  (44)  of  I]?|_(/i),  where  p®  =  0  and  the 
trajectory  ‘jumps’  to  S(_(/i)  (the  right  panels  in  Fig.  8)  after  which  it  returns  back  deeper  into  the  reactants 
region  with  px  <  0. 

Figure  9  shows  the  image  of  a  Lagrangian  cylinder  Aj  j  for  /  <  0  together  with  a  trajectory  contained 
in  Aj!"j.  The  trajectory  approaches  the  bottleneck  region  with  negative  Px  from  the  products  side  of  the 
energy  surface  within  SL(h)  (the  right  panels  of  Fig.  9),  reaches  the  boundary  (44)  of  SL(h),  where  Px  =  ^ 
and  the  trajectory  ‘jumps’  to  S(|_(/i)  (the  left  panels  in  Fig.  9),  in  which  it  returns  back  deeper  in  to  the 
products  region  with  >  0. 

Figure  10  shows  the  F^  images  of  Lagrangian  cylinders  Aj  j  and  Aji'j  for  /  >  0.  For  /  >  0,  Aj  j  contains 
forward  reactive  trajectories  which  have  Px  >  0  and  therefore  are  completely  contained  in  Sj|_(/i)  (the  left 
panels  in  Fig.  10),  and  Aj j  contains  backward  reactive  trajectories  which  have  Px  <  0  and  therefore  are 
completely  contained  in  SL(h)  (the  right  panels  in  Fig.  10).  The  forward  and  backward  trajectories  spiral 
about  the  forward  and  backward  reaction  paths  which  are  also  shown  in  Fig.  10. 

The  images  of  the  dividing  surface,  Sjg(/i),  and  the  NHIM,  Snhim(^);  under  F^  are  shown  in  Fig.  11. 
The  forward  hemisphere  f{h)  is  contained  in  S5|_(/i)  (the  left  panels  in  Fig.  11);  the  backward  hemisphere 
^ds  b(^)  contained  in  SL(h)  (the  right  panels  in  Fig.  11).  The  periodic  orbit  §nhim(^)  which  forms  the 
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Figure  10:  Model  I:  forward  and  backward  reactive  trajectories  spiraling  about  the  forward  and  backward 
reaction  paths,  respectively.  The  yellow  wireframe  in  the  top  left  panel  indicates  a  Lagrangian  cylinder 
with  /  >  0  contained  in  E5|_(/i);  the  blue  wireframe  in  the  top  right  panel  indicates  a  Lagrangian  cylinder 
Aj  j  with  /  >  0  contained  in  SL  (/*-)■  The  black  lines  show  trajectories  contained  in  the  respective  cylinders. 
The  blue  centerlines  show  the  forward  (left)  and  backward  (right)  reaction  paths.  The  bottom  panels  show 
the  intersections  of  the  Lagrangian  cylinders  with  the  plane  y  =  0  together  with  the  forward  and  backward 
reaction  paths. 
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Figure  11:  Model  I:  the  dividing  surface  and  the  NHIM.  The  dark  red  wireframes  indicate  the  two  hemispheres 
f{h)  (top  left  panel)  and  ^^{h)  (top  right  panel)  the  dividing  surface  The  light  blue  circle  in  the 

top  panels  shows  the  periodic  orbit  Snhim(^)-  The  lower  panels  show  the  intersections  of  f{h),  B'^^  ^^{h) 
and  §nhim(^)  with  the  plane  y  =  0- 
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Figure  12:  Model  I:  the  stable  and  unstable  manifolds  of  the  NHIM.  The  wireframes  in  the  top  left  panel 
show  the  branches  IF/  (green)  and  W/  (red),  and  the  wireframes  in  the  top  right  panel  show  the  branches 
IF/  (green)  and  W/*  (red)  of  the  stable  and  unstable  manifolds  of  Sj,jjjj^(/i).  For  each  branch,  a  trajectory 
(black  line)  contained  in  the  respective  branch  is  shown.  The  bottom  panels  show  the  intersections  of  the 
branches  with  the  plane  y  =  0. 


equator  of  Sds(^)  Px  =  0  and  thus  runs  along  the  boundaries  of  Sj|_(/i)  and  EL(/i). 

Figure  12  shows  the  images  of  the  stable  and  unstable  manifolds  VF®  and  1F“  of  Snhim(^)-  "^he  forward 
branches  W/  and  IF/  are  contained  in  S(|_(/i)  (left  panels  in  Fig.  12),  and  the  backward  branches  IF/  and 
IF/  are  contained  in  SL(/i)  (left  panels  in  Fig.  12).  Figure  12  also  shows  how  trajectories  contained  in  the 
branches  asymptotically  approach  the  periodic  orbit  §nhim(^)  respective  time  limit  t  ±oo. 

Regarding  Fig.  12  together  with  Fig.  10  we  see  that  the  forward  reactive  cylinder  W/  U  IF/  contained  in 
E5|_(/i)  encloses  all  forward  reactive  trajectories,  and  the  backward  reactive  cylinder  IF/U  W/  encloses  all  the 
backward  reactive  trajectories.  Note  that  the  images  of  the  forward  and  backward  reactive  cylinders  are 
smooth  along  their  junctions  given  by  the  periodic  orbit  Si;jhim(^)  while  in  the  original  phase  space  IF/  U  IF/ 
and  IF/  U  IF/  have  a  kink  along  §i;jhim(^)'  This  is  an  artifact  of  the  map  F^  which,  as  mentioned  above,  is 
singular  along  the  cylinder  Ua;g]RSj  ^{h)  which  contains  Sj;jjjj]yj(/i)  =  Sj  o(^)- 

3.2  Model  II:  stereographic  projection  of  2-spheres 

The  main  idea  of  the  second  model  (Model  II)  is  to  represent  the  2-spheres  S^(/i)  defined  by  Eq.  (23)  in  terms 
of  stereographic  type  projections  to  a  plane.  Recall  that  §^(/i)  is  given  by  the  axisymmetric  ellipsoid  (23). 
We  define  a  stereographic  projection  of  this  ellipsoid  by  projecting  each  point  on  it  along  the  line  through 
this  point  and  the  South  pole  (px^Py^y)  =  (— -\/(2/i  +  Xx‘^)/X,  0, 0)  of  the  ellipsoid  to  the  equatorial  plane 
of  the  ellipsoid.  This  projection  is  illustrated  in  Fig.  13.  Points  on  the  hemisphere  ^{h)  of  S^(/i),  which 
we  defined  in  (37),  are  projected  to  the  circular  disk  bounded  by  the  equator  (39)  (e.g.,  the  point  Pi  in 
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Figure  13:  Model  II:  stereographic  projection  of  to  its  equatorial  plane  with  coordinates  {y,py).  The 
points  Pi  and  P2  are  the  images  of  the  points  Pi  and  P2  under  the  projection  (47). 


Fig.  13).  The  North  pole  is  mapped  to  the  center  of  this  disk.  Points  on  the  hemisphere  ^{h)  defined  in 
(38)  are  projected  to  points  in  the  plane  outside  of  the  disk  (see,  e.g.,  the  point  P2  in  Fig.  13).  The  South 
pole  {pxiPy,y)  =  (— \/(2/i  +  Ax2)/A,  0, 0)  is  mapped  to  infinity.  The  equator  S()  ,j,(/i)  of  §^(/i)  (see  (39))  is 
invariant  under  the  projection.  Formally  this  projection  is  given  by 

TTsgx  :  S^(/l)  ^  IR.^  {Px,Py,y)  ^  {y,Py)  ,  (47) 

where 


_ 2 

and  E,  denotes  the  compactified  plane  U  {oo}.  Note  that  the  familiar  standard  stereographic  projection 
is  from  a  2-sphere  of  unit  radius  to  its  equatorial  plane  where  the  ‘base  point’  of  the  projection  is  chosen  to 
be  the  North  pole  rather  than  the  South  pole  like  in  our  choice.  This  explains  the  scaling  factors  and  signs 
in  (48)  which  differ  from  the  standard  formulae.  We  use  the  projection  TTggx  in  (47)  to  define  a  map  of  points 
(x,y,px,Py)  on  the  energy  surface  S(/i)  according  to 

-F”  :  {x,y,Px,Py)  {x,7rsgxiPx,Py,y))  =  {x,y,Py).  (49) 

The  map  is  singular  along  the  North  poles  of  the  spheres  S^(/i),  a;  G  E,  where  y  and  py  are  zero,  and 
Px  reaches  its  minimal  value  for  fixed  x  and  h.  The  set  of  singular  points  thus  is  the  backward  reaction 
path  which  is  mapped  to  infinity.  As  a  result  the  image  of  S(/i)  under  is  the  Cartesian  product  E  x  E  . 
Note  that  the  map  preserves  the  symmetry  (34),  or  more  precisely,  like  in  Model  I,  manifolds  which  are 
invariant  under  the  symmetry  action  (34)  are  rotationally  symmetric  about  the  x  axis  of  the  image  of 

In  the  following  we  use  this  representation  to  show  geometrical  visualizations  of  the  various  phase  space 

_ 2 

structures  that  govern  reaction  dynamics  analogously  to  Sec.  3.1.  Since  the  energy  surface  maps  to  E  x  E 
under  F^^  we  omit  a  picture  of  the  image  of  the  energy  surface.  In  Fig.  14  we  show  the  F^^  image  of  a 
Lagrangian  manifold  Aj  j  with  /  <  0  together  with  a  nonreactive  trajectory  contained  in  this  cylinder.  The 
trajectory  comes  in  from  the  reactants  component  (x  <C  —1)  on  the  left  with  positive  Px,  spiraling  about  the 
X  axis.  The  radius  of  the  rotation  increases  as  x  increases.  After  reaching  a  maximal  x  value  (where  Px  =  0) 
the  trajectory  returns  deeper  and  deeper  into  the  reactants  region,  still  spiraling  about  the  x  axis  with  the 
radius  of  rotation  increasing  as  x  decreases.  In  Fig.  15  we  show  the  analogous  picture  for  the  F^^  image 
of  a  Lagrangian  manifold  A)!'j  with  /  <  0.  The  trajectories  come  in  from  the  products  component  on  the 
right  (x  ^  1)  with  negative  Px  undergoing  rotations  with  increasing  radius  about  the  x  axis  as  x  decreases. 


Pyy/  {2h  +  Xx‘^)/X 
^  (2h  +  Xx"^)  /  X  +  Px 


j/a/  {2h  +  Xx^)/X 
/{2h  +  Xx'^)/X  +Px 
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Figure  14:  Model  II:  nonreactive  trajectories  that  stay  in  the  reactants  component.  The  orientation  of  the 
X,  y  and  Py  axes  is  indicated  by  the  dreibein.  The  wireframe  in  the  left  panel  shows  a  Lagrangian  cylinder 
Ajj  defined  in  (14)  for  /  <  0.  The  black  line  shows  one  of  the  trajectories  contained  in  Aj  j.  The  red  lines 
in  the  right  panel  show  the  intersection  of  Aj  j  with  the  plane  y  =  0. 


Figure  15:  Model  II:  nonreactive  trajectories  that  stay  in  the  products  component.  The  green  wireframe  in 
the  left  panel  shows  a  Lagrangian  cylinder  A~^  j  defined  in  (13)  for  I  <  0.  The  black  line  shows  one  of  the 
trajectories  contained  in  Af  j.  The  green  lines  in  the  right  panel  show  the  intersection  of  Af  j  with  the  plane 
y  =  0. 
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Figure  16:  Model  II:  forward  reactive  trajectories  spiraling  about  the  forward  reaction  path.  The  yellow 
wireframe  in  the  left  panel  indicates  a  Lagrangian  cylinder  Aj!"  j  with  /  >  0.  The  black  line  shows  a  trajectory 
contained  in  Aj  j.  The  blue  line  shows  the  forward  reaction  path.  The  right  panel  shows  the  intersection  of 
Af  j  with  the  plane  y  =  0  together  with  the  forward  reaction  path. 


Figure  17:  Model  II:  backward  reactive  trajectories.  The  blue  wireframe  in  the  left  panel  indicates  a  La¬ 
grangian  cylinder  Ajj  with  /  >  0.  The  black  line  shows  a  trajectory  contained  in  A^j.  The  right  panel 
shows  the  intersection  of  Ay  j  with  the  plane  y  =  0. 


After  reaching  a  minimal  x  value,  the  trajectory  returns  deeper  and  deeper  into  the  products  region,  still 
undergoing  rotation  about  the  x  axis,  where  the  radius  of  the  rotation  now  increases  as  x  increases. 

Figure  16  shows  the  image  of  a  Lagrangian  cylinder  Aj  j  with  /  >  0  together  with  the  forward 
reaction  path  which  coincides  with  the  x  axis.  Trajectories  contained  in  A'f  j  evolve  from  the  reactants  to 
the  products  spiraling  about  the  forward  reaction  path.  Figure  17  shows  the  image  of  a  Lagrangian 
cylinder  Ay j  with  J  >  0  that  contains  backward  reactive  trajectories.  These  trajectories  spiral  with  much 
larger  radii  of  rotation  than  the  forward  reactive  trajectories  in  Fig.  16.  In  fact,  since  the  backward  reaction 
path  is  mapped  to  the  line  R  U  cx)  the  backward  reactive  trajectories  can  again  be  considered  to  spiral  about 
the  backward  reaction  path. 

In  Fig.  18  we  show  the  images  of  the  dividing  surface  and  the  NHIM  §nhim(^)  under  The 

dividing  surface  maps  to  the  compactified  plane  a;  =  0.  The  F^^  image  of  the  NHIM,  §nhim(^)i  divides 
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Figure  18:  Model  II:  dividing  surface  and  NHIM.  The  dark  red  wireframe  in  the  left  panel  shows  the  dividing 
surface  The  periodic  orbit  Snhim(^)  (blue  circle)  divides  it  into  the  forward  hemisphere  ^{h) 

(inner  circular  disk)  and  the  backward  hemisphere  ^{h)  (outer  infinite  annulus).  The  right  panel  shows 
the  intersection  of  Sds(ft.)  and  Snhim(^)  the  plane  y  =  0- 


Figure  19:  Model  II:  forward  reactive  cylinder,  forward  dividing  surface  and  NHIM.  The  green  and  red 
wireframes  in  the  left  panel  indicate  the  branches  Wf  and  IF"  of  the  stable  resp.  unstable  manifolds  of  the 
NHIM  Si;jhim(^)  (blue  circle).  The  black  lines  show  trajectories  contained  in  Wf  and  .  The  dark  red 
wireframe  at  the  center  which  is  difficult  to  see  shows  the  forward  hemisphere  B^^  ^(h)  of  the  dividing  surface. 
The  right  panel  shows  the  intersections  of  the  structures  on  the  right  with  the  plane  y  =  0. 

it  into  the  forward  hemisphere  B^^  f{h)  which  appears  as  a  circular  disk  in  the  plane  a;  =  0  centered  about 
the  origin  and  the  backward  hemisphere  B^^  hW  "'^hich  appears  as  an  annulus  that  extends  to  the  boundary 
of  at  infinity.  The  image  of  the  forward  reactive  trajectories  in  Fig.  16  which  are  located  near  the  x 
axis  cross  the  image  of  B^^  f{h),  and  the  F^^  image  of  the  backward  reactive  trajectories  in  Fig.  17  which 
stay  away  from  the  x  axis  cross  the  F^^  image  of  b(^)'  The  two  types  of  trajectories  are  enclosed  by  the 
forward  and  backward  reactive  cylinders  U  Wf“  and  U  whose  F^^  images  are  shown  in  Figs.  19 
and  20,  respectively. 

Note  that  in  this  geometrical  representation  of  the  energy  surface  (Model  H)  the  bottleneck  area  of  S(/i) 
associated  with  the  reaction  region  containing  the  dividing  surface  is  not  so  pronounced  as  in  Model  I. 
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Figure  20:  Model  II:  backward  reactive  cylinder,  backward  dividing  surface  and  NHIM.  The  green  and  red 
wireframes  in  the  left  panel  indicate  the  branches  W-^  and  of  the  stable  resp.  unstable  manifolds  of 
the  NHIM  §nhim(^)  (blue  circle).  The  black  lines  show  trajectories  contained  in  and  W^.  The  dark 
red  wireframe  shows  the  backward  hemisphere  ^{h)  of  the  dividing  surface.  The  right  panel  shows  the 
intersections  of  the  structures  on  the  right  with  the  plane  y  =  0- 


3.3  Model  III:  McGehee  representation 

The  main  idea  of  the  third  model  (Model  III)  of  the  energy  surface  S(h)  is  to  represent  the  family  of 
(topological)  2-spheres,  §^(/i),  in  Eq.  (23)  parametrized  by  x  €  IR  as  a  ‘nested’  set  of  concentric  2-spheres  in 
R^.  This  can  be  achieved  by  projecting  each  of  the  2-spheres,  E>‘^{h)  to  a  2-sphere  of  a  radius  which  depends 
on  X.  One  way  to  realize  this  is  to  map  a  point  P  =  {px,Py,  y)  on  the  §^(h)  according  to 

s  (^) 

TTMcGx  :  {y,Px,Py)  ^  iy,Px,Py)Y( - pi  =’  iy^Px,Pv)  ■  (50) 

II  [yjPxyPyJW 

Here  H-H  denotes  the  Euclidean  norm,  i.e.  \\{px,Py,y)\\  =  (p^  +  Py  +  and  s  is  the  scaling  function 

s(x)  :=  1 -I- i  tanhx  (51) 

which  maps  the  domain  x  G  R  monotonicly  to  the  interval  (1/2,  3/2).  The  projection  ttmcGx  is  illustrated  in 
Fig.  21.  We  name  this  projection  after  McGehee  who  first  used  this  type  of  construction  to  present  an  energy 
surface  of  the  structure  x  R  [McGehee(1969)].  A  similar  construction  can  also  be  found  in  the  work  by 
MacKay  [MacKay(1990)]. 

We  use  ttmcGx  in  (50)  to  define  a  map  of  points  {x,y,Px,Py)  on  the  energy  surface  E(h)  to  R^ 
according  to 

:  {x,y,Px,Py)  ^  TTMcGxiPx,  Py,y)  =■  {y,Px,Py)  ■  (52) 

Like  the  maps  and  the  map  also  preserves  the  symmetry.  In  Model  HI,  manifolds  that  are 
invariant  under  the  symmetry  action  (34)  have  F^^^  images  that  are  symmetric  about  the  Px  axis. 

The  image  of  the  energy  surface  E(/i)  of  topology  R  x  under  the  map  F^^^  is  the  ‘spherical  shell’ 
consisting  of  the  open  ball  of  radius  3/2  in  R^  minus  the  closed  ball  of  radius  1/2.  This  image  is  shown  in 
Fig.  22.  The  reactants  region  x  <C  —1  is  contained  in  the  outer  region  of  this  spherical  shell;  the  products 
region  x  »  1  is  contained  in  the  inner  region.  In  Fig.  22  we  also  marked  the  points  {y,Py,Px)  =  (0,0,  ±1/2) 
and  {y,Py,Px)  =  (0,  0,  ±3/2)  to  which  the  F^^^  images  of  all  trajectories  converge  as  time  goes  to  — oo  or  ±oo 
(as  long  as  the  trajectory  is  unbounded  in  the  respective  limit).  This  is  a  consequence  of  the  conservation  of 
the  integrals  I  and  J  (the  energy  in  the  x  and  y  degrees  of  freedom,  respectively)  which  implies  that  y  and 
Py  are  bounded  for  each  trajectory  whereas  Px  diverges  as  x  goes  to  — oo  or  ±oo.  In  Fig.  23  we  show  the 
example  of  the  F^^^  image  of  a  Lagrangian  cylinder  Aj j  defined  for  /  <  0.  The  trajectories  in  Ajj  come  in 
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Figure  21:  Model  III:  McGehee  projection  of  a  point  P  =  {y,Px,Py)  on  the  (topological)  2-sphere  S^(/i) 
(inner  dark  shaded  sphere)  to  the  point  P  =  P  =  {y,Px,Py)  on  the  sphere  of  radius  s{x)  (outer  light  shaded 
sphere)  with  s{x)  defined  in  (51).  Note  that,  depending  on  x,  the  outer  sphere  could  also  be  inside  of  the 
inner  sphere,  or  the  two  could  overlap. 


from  deep  in  the  reactants  region  where  x  <IC  —1  and  Px  3>  1  which  is  asymptotically  (in  the  limit  t  —>■  — oo) 
mapped  to  {y,Py,Px)  =  (0,0,  3/2)  under  and  return  deep  into  the  reactants  region  where  a;  <C  — 1  and 
Pa;  <C  —  1  which  is  asymptotically  (in  the  limit  t  oo)  mapped  to  {y,Py,Px)  =  (0,  0,  —3/2)  under 

In  Fig.  24  we  show  the  F^^^  image  of  a  Lagrangian  cylinder  A/j  with  I  <  0  which  contains  nonreactive 
trajectories  in  the  products  region.  They  come  in  from  deep  in  the  products  region  x  1  with  pa,  <C  —  1 
which  under  F^^^  is  mapped  to  {y,Py,Px)  =  (0, 0,  —1/2)  and  move  back  deep  into  the  products  region  x  1 
with  Pa;  1  which  under  F^^^  is  mapped  to  {y,Py,Px)  =  (0,  0, 1/2). 

Fig.  25  shows  the  F^^^  image  of  a  Lagrangian  cylinder  A/j  which  has  /  >  0.  The  Lagrangian  manifold 
consists  of  forward  reactive  trajectories  which  spiral  about  the  forward  reaction  path  that  is  also  shown  in 
Fig.  25.  Fig.  26  shows  the  analogous  image  of  a  Lagrangian  cylinder  Aj  j  which  has  /  >  0  and  hence  consists 
of  backward  reactive  trajectories. 

For  clarity  we  show  the  F^^^  images  of  the  branches  IF/,  IF/,  IF/  and  IF/  of  the  stable  and  unstable 
manifolds  of  the  NHIM  separately  in  Figs.  28-31.  The  F^^^  image  of  the  union  IF/  U  IF/  forms  the  forward 
reactive  cylinder  in  Model  III.  It  encloses  the  F^^^  image  of  all  forward  reactive  trajectories  as  can  be  seen 
from  comparing  Fig.  25  with  Fig.  32.  Similarly  the  F^^^  image  of  the  backward  reactive  cylinder  IF/  U  IF/ 
encloses  the  F^^^  image  of  all  backward  reactive  cylinder  as  can  be  seen  from  comparing  Fig.  26  with  Fig.  33. 

3.4  Implications  for  Nonlinear  Hamiltonian  Vector  Fields 

We  have  chosen  to  develop  the  three  models  of  the  phase  space  structures  near  a  saddle  x  centre  equilibrium  by 
considering  a  quadratic,  two  degree  of  freedom  Hamiltonian  which  is  the  normal  form,  through  second  order, 
in  the  neighborhood  of  such  an  equilibrium  point.  The  advantage  of  this  is  that  is  provided  very  explicit 
formulae  for  the  phase  space  structures  that  allowed  straightforward  graphical  representation.  However, 
precisely  the  same  geometrical  structures  exist  for  nonlinear  Hamiltonian  vector  fields  in  the  neighborhood  of 
a  saddle  x  center  equilibrium  point.  In  the  nonlinear  case  the  phase  space  structures  can  again  be  described 
by  explicit  mathematical  formulae  using  a  normal  expansion  at  the  saddle  x  center  equilibrium  (see  Sec.  2.2, 
and  [Waalkens  et  al.(2008)]  for  a  detailed  discussion).  For  a  given  system,  this  in  general  leads  a  Hamilton 
function  which  is  a  power  series  rather  than  a  linear  function  of  the  integrals  I  and  J  (see  (20)).  The 
coefficients  of  this  powers  series  are  obtained  from  the  normal  form  calculations  for  an  explicit  systems. 

Nevertheless,  the  normal  form  theory  allows  us  to  conclude  that  the  same  geometrical  structures,  having 
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Figure  22:  Model  III:  the  energy  surface.  The  left  panel  shows  the  energy  surface  in  the  threedimensional 
space  with  coordinates  (y,Py,Px)-  The  orientation  of  the  y,  Py  and  px  axes  are  indicated  by  the  dreibein. 
(Note  that  the  center  of  the  coordinate  system  is  not  the  center  of  the  dreibein.)  The  black  points  mark  the 
points  {y,Py,Px)  =  (0,  0,  ±1/2)  and  {y,Py,Px)  =  (0, 0,  ±3/2)  (see  text).  The  right  panel  shows  the  intersection 
of  the  energy  surface  with  the  plane  y  =  0. 


Figure  23:  Modell  III:  nonreactive  trajectories  that  stay  in  the  reactants  component.  The  red  wireframe  in 
the  left  panel  indicates  a  Lagrangian  cylinder  Aj  j  defined  in  (14)  for  I  <  0.  The  black  line  shows  one  of  the 
trajectories  contained  in  Aj  j.  The  red  lines  in  the  right  panel  show  the  intersection  of  Aj  j  with  the  plane 
y  =  0. 
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Figure  24:  Modell  III:  nonreactive  trajectories  that  stay  in  the  products  component.  The  green  wireframe 
in  the  left  panel  indicates  a  Lagrangian  cylinder  A'^  j  defined  in  (13)  for  /  <  0.  The  black  line  shows  one  of 
the  trajectories  contained  in  A'^j.  The  green  lines  in  the  right  panel  show  the  intersection  of  A~^  j  with  the 
plane  y  =  0. 


Figure  25:  Model  III:  forward  reactive  trajectories  spiraling  about  the  forward  reaction  path.  The  yellow 
wireframe  in  the  left  panel  indicates  a  Lagrangian  cylinder  A'j^  j  with  /  >  0.  The  black  line  shows  a  trajectory 
contained  in  Aj  j.  The  blue  line  shows  the  forward  reaction  path.  The  right  panel  shows  the  intersection  of 
Aj!"  j  with  the  plane  y  =  0  together  with  the  forward  reaction  path. 
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Figure  26:  Model  III:  backward  reactive  trajectories  spiraling  about  the  backward  reaction  path.  The  blue 
wireframe  in  the  left  panel  indicates  a  Lagrangian  cylinder  Aj"  j  with  /  >  0.  The  black  line  shows  a  trajectory 
contained  in  Aj" j  spiraling  about  the  backward  reaction  path  marked  in  blue.  The  right  panel  shows  the 
intersection  of  Aj"  j  with  the  plane  y  =  0  together  with  the  backward  reaction  path. 


Figure  27:  Model  III:  dividing  surface  and  NHIM.  The  dark  red  wireframe  in  the  left  panel  shows  the  dividing 
surface  S‘^^{h).  The  periodic  orbit  Snhim(^)  (blue  circle)  divides  it  into  the  forward  hemisphere  f{h)  (the 
hemisphere  on  the  right)  and  the  backward  hemisphere  ^^{h)  (the  hemisphere  on  the  left).  The  right 
panel  shows  the  intersection  of  and  §nhim(^)  '"^bh  the  plane  y  =  0- 
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Figure  28:  Model  III:  the  forward  branch  of  the  stable  manifold.  The  green  wireframe  in  the  left  panel 
shows  Wf  .  The  black  line  shows  a  trajectory  contained  in  W^.  The  green  lines  in  the  right  panels  show  the 
intersection  of  with  the  plane  y  =  0- 


Figure  29:  Model  III:  the  forward  branch  of  the  unstable  manifold.  The  red  wireframe  in  the  left  panel 
shows  IFf“.  The  black  line  shows  a  trajectory  contained  in  Wf“.  The  red  lines  in  the  right  panels  show  the 
intersection  of  Wf*  with  the  plane  y  =  0- 
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Figure  30:  Model  III:  the  backward  branch  of  the  stable  manifold.  The  green  wireframe  in  the  left  panel 
shows  W^.  The  black  line  shows  a  trajectory  contained  in  W^.  The  green  lines  in  the  right  panels  show  the 
intersection  of  with  the  plane  y  =  0. 


Figure  31:  Model  III:  the  backward  branch  of  the  unstable  manifold.  The  red  wireframe  in  the  left  panel 
shows  W^.  The  black  line  shows  a  trajectory  contained  in  W^.  The  red  lines  in  the  right  panels  show  the 
intersection  of  with  the  plane  ji  =  0. 
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Figure  32:  Model  III:  the  forward  reactive  cylinder,  the  dividing  surface  and  the  NHIM  as  a  superposition  of 
Figs.  27,  28  and  29. 


Figure  33:  Model  III:  the  backward  reactive  cylinder,  the  dividing  surface  and  the  NHIM  as  a  superposition 
of  Figs.  27,  30  and  31. 
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the  same  interpretation  and  meaning,  hold  for  more  general,  nonlinear  Hamiltonian  vector  fields  in  the 
neighborhood  of  a  saddle  center-. .  .-center  equilibrium  point,  but  obtaining  the  explicit  formulae  will  require 
a  normal  form  calculation  for  a  given  system.  It  is  in  this  sense  that  restricting  ourselves  to  quadratic 
Hamiltonians  is,  indeed,  ’’without  loss  of  generality”,  and  it  provides  the  method  and  approach  for  visualizing 
the  same  phase  space  structures  governing  reaction  dynamics  in  more  general  systems. 


4  A  model  of  the  n-DoF  system  in  the  space  of  the  integrals  of 
motion 

For  a  system  with  more  than  2  degrees  of  freedom,  the  dimensionality  of  the  phase  space  structures  discussed 
in  Sec.  2  is  too  high  to  allow  for  a  similar  explicit  visualization  like  for  the  2-DoFmodels  in  Sec.  3.  Instead 
of  projecting  the  phase  space  structures  to  the  planes  of  the  coordinates  and  conjugate  momenta  of  the 
normal  form  like  in  Sec.  2.2  it  is  useful  to  present  the  phase  space  structures  for  n  DoFin  the  space  of  the 
integrals  I,  J2,  ■  ■  ■ ,  Jn-  This  can  be  accomplished  using  the  momentum  map  A4  defined  in  (12)  in  Sec.  2.3. 
The  discussion  here  follows  [Waalkens  et  al.(2008)]. 

As  mentioned  in  Sec.  2.3  a  fibre  (the  joint  level  set  of  the  integrals  /,  J2,...,  Jn  in  phase  space)  is 
called  singular  if  it  contains  one  (or  more)  irregular  point(s).  The  image  of  all  the  singular  fibres  under  the 
momentum  map  is  called  the  bifurcation  diagram.  It  is  easy  to  see  that  the  bifurcation  diagram  consists  of 
the  set  of  (/,  J2, . . . ,  Jn),  where  one  or  more  of  the  integrals  vanish.  In  Fig.  34  we  show  the  image  of  the 
energy  surface  with  energy  h  >  0  under  the  momentum  map  M  in  the  space  of  the  integrals  for  n  =  3  degrees 
of  freedom. 


Figure  34:  Sketch  of  the  image  of  the  energy  surface  of  energy  h  >  0  under  the  momentum  map  M  in 
Equation  (12)  in  the  space  of  the  integrals  I  and  Jk,  k  =  2, . . .  ,n,  for  the  case  of  n  =  3  degrees  of  freedom. 
The  green/dark  pink  brindled  piece  of  the  image  of  the  energy  surface  has  /  >  0;  the  turquoise  piece  has 
/  <  0.  The  intersections  with  the  planes  /  =  0,  J2  =  0  and  J3  =  0  (pieces  of  which  are  visualised  by 
semitransparent  planes  for  clarity)  form  the  bifurcation  diagram  of  the  energy  surface.  The  image  of  the 
energy  surface  is  not  bounded  in  the  direction  of  negative  I  as  indicated  by  the  dashed  line  at  the  bottom. 
The  topology  of  the  fibres  J2,  J3)  is  indicated  for  the  various  points  (/,  J2,  J3)  marked  by  dots.  The 

fibre  of  a  point  (7,  J2,  J3)  with  7  7^  0  consists  of  two  disconnected  manifolds  as  indicated  by  the  factor  of  2. 
The  fibre  of  a  point  (7,  J2,  J3)  with  7  =  0  consists  of  a  single  connected  manifold. 


The  bifurcation  diagram  (of  the  energy  surface)  consists  of  the  intersections  of  the  image  of  the  energy 
surface  (the  turquoise  and  green/dark  ping  brindled  surface  in  Fig.  34)  with  one  of  the  planes  7  =  0,  J2  =  0 
or  J3  =  0.  Upon  approaching  one  of  the  edges  that  have  J2  =  0  or  J3  =  0  the  circle  in  the  plane  (92, 7*2)  or 
(93, Pa),  respectively,  shrinks  to  a  point,  and  accordingly  the  regular  fibres  x  H  which  form  the  Lagrange 
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manifolds  and  (see  (13)  and  (14))  reduce  to  cylinders  or  ‘tubes’  x  R.  At  the  top 

corner  in  Fig.  34  both  J2  and  J3  are  zero.  Here  both  circles  in  the  centre  planes  (92, P2)  and  {q^iPz)  have 
shrinked  to  points.  The  corresponding  singular  fibre  consists  of  two  lines,  R,  which  are  the  forward  and 
backward  reaction  paths,  respectively  (see  also  Fig.  4). 

The  fibres  mentioned  so  far  all  have  /  7^  0  and  each  consist  of  a  pair  of  two  disconnected  components. 
For  /  <  0,  one  member  of  each  pair  is  located  on  the  reactants  side  and  the  other  on  the  products  side  of 
the  dividing  surface.  For  /  >  0,  one  member  of  each  pair  consists  of  trajectories  evolving  from  reactants  to 
products  and  the  other  member  consists  of  trajectories  that  evolve  from  products  to  reactants.  In  fact  the 
two  members  of  a  fibre  which  has  /  >  0  are  contained  in  the  energy  surface  volume  enclosed  by  the  forward 
and  backward  reactive  spherical  cylinders  Wf{h)  and  H4(/i),  see  Fig.  4.  For  this  reason  we  marked  the  piece 
of  the  image  of  the  energy  surface  under  the  momentum  map  which  has  /  >  0  by  the  same  green/dark  pink 
colour  in  Fig.  34  that  we  used  Fig.  4.  Green  corresponds  to  forward  reactive  trajectories  and  dark  pink 
corresponds  to  backward  reactive  trajectories.  Under  the  momentum  map  these  trajectories  have  the  same 
image. 

The  light  blue  line  in  Fig.  34  which  has  /  =  0  is  the  image  of  the  NHIM  under  the  momentum  map.  For 
three  degrees  of  freedom,  the  NHIM  is  a  3-dimensional  sphere,  and  as  we  will  discuss  in  more  detail  in  Sec.  5 
and  indicated  in  Fig.  34  it  is  foliated  by  a  one-parameter  family  of  invariant  2-tori  which  shrink  to  periodic 
orbits,  i.e.  circles,  S^,  at  the  end  points  of  the  parameterisation  interval. 


5  Structure  of  the  NHIM  for  n-DoF 

We  have  described  how  the  NHIM,  its  stable  and  unstable  manifolds,  and  the  dividing  surface,  sit  within  the 
phase  space  and  govern  the  reaction  dynamics  in  several  2DoF  geometrical  models.  However,  we  have  not 
given  any  indication  of  the  structure  of  the  NHIM  itself,  except  in  the  2-DoF  case,  in  which  the  NHIM  is 
a  single  periodic  orbit.  In  this  section  we  discuss  various  ways  of  visualizing  the  geometrical  structure  and 
dynamics  associated  with  the  NHIM. 

Recall  from  Sec.  2.2  that  in  terms  of  the  local  decoupling  provided  by  the  normal  form  into  reaction  coor¬ 
dinates  and  bath  modes,  points  on  the  NHIM,  SNfflM(^)>  contain  zero  ‘energy’  in  the  reactive  (saddle)  mode 
(i.e.,  1  =  0);  all  ‘energy’  is  contained  in  the  bath  (centre)  modes,  which  results  in  oscillatory,  quasiperiodic 
motion  with  n  —  1  independent  frequencies.  As  also  mentioned  in  Sec.  2.2,  this  leads  to  the  interpretation  of 
the  NHIM  as  the  energy  surface  of  an  invariant  subsystem  with  one  degree  of  freedom  less  than  the  full  system 
(a  supermolecule  located  for  a  frozen  reaction  coordinate  (and  conjugate  momentum)  between  reactants  and 
products). 

Since,  typically  (for  more  than  1  DoF),  there  is  more  than  one  bath  mode  present,  the  total  energy 
can  be  distributed  in  a  variety  of  ways  amongst  the  bath  modes.  For  example,  putting  all  of  the  energy 
into  any  one  of  the  bath  modes  results  in  a  single  periodic  orbit  (often  referred  to  as  a  Lyapunov  orbit). 
More  generally,  the  supermolecule  of  fixed  energy  h  has  energy  distributed  between  all  of  the  n  —  1  bath 
modes.  Since  the  integrals  of  motion  are  constant  on  trajectories,  the  amount  of  energy  in  each  particular 
bath  mode  remains  fixed  during  the  motion.  Such  a  molecule  therefore  undergoes  quasiperiodic  motion;  it 
oscillates  independently  in  the  n  —  1  different  bath  modes.  The  motion  of  a  typical  trajectory  in  the  NHIM  is 
quasiperiodic  on  an  (n  —  l)-torus.  Heuristically,  one  can  think  of  the  NHIM  as  the  invariant  surface  made  up 
of  trajectories  with  all  possible  energy  distributions  between  the  bath  modes  (and  no  energy  in  the  reactive 
mode). 

5.1  The  NHIM  („  =  3  DoF) 

As  noted  above,  the  NHIM  for  2-DoF  systems  is  simply  a  periodic  orbit.  This  is  an  immediate  consequence 
of  the  fact  that  there  is  only  a  single  bath  (or  centre)  mode.  In  this  section  we  consider  the  next  most 
complex  case  in  some  detail,  the  3-DoF  situation  (2  bath  modes),  for  which  the  NHIM  is  diffeomorphic  to 
S^.  However,  first  we  will  describe  a  way  in  which  can  be  visualized  in  lower  dimensions.  This  is  helpful 
from  the  point  of  view  of  visual  intuition.  Our  discussion  in  this  section  owes  a  great  deal  to  the  nice  report 
of  [Chisholm(2000)]. 
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Geometrical  Construction  of  in  Three  Dimensions.  To  begin  with,  first  note  that  in  the  usual 
manner,  to  draw  a  1-sphere  (i.e.  a  circle),  we  need  two  dimensons,  to  draw  a  2-sphere  we  need  three 
dimensions,  and  to  draw  a  3-sphere  we  need  four  dimensions.  We  will  describe  a  method  of  drawing  the 
3-sphere  which  requires  only  three  dimensions.  This  particular  representation  will  enable  us  to  make  a  clear 
picture  of  the  dynamics  on  the  NHIM.  The  method  can  be  viewed  as  a  generalization  of  the  presentation 
of  a  two-dimensional  sphere,  S^,  in  terms  of  two  hemispheres  (topological  2-balls,  B^)  and  an  equator  (a 
onedimensional  sphere,  S^)  that  we  discussed  in  Sec.  3.1.  For  S^,  one  can  form  a  similar  representation  by 
using  two  open  3-balls,  to  represent  “Northern”  and  “Southern”  hemispheres,  and  identifying  points  on 
their  boundaries,  S^,  to  make  an  “equator”; 

=  i?North  U  §|quator  ^  .Bsouth’  (53) 

We  will  choose  the  “North”  pole  (N)  to  be  at  the  centre  of  the  North  ball  and  the  “South”  (S)  pole  at  the 
centre  of  the  South  ball.  This  model  is  illustrated  in  Fig.  35. 


Figure  35:  A  model  of  the  sphere  formed  by  taking  two  open  3-balls,  i^North  -^South’  identifying 
points  on  their  boundaries;  the  “equator”  Slq^ator-  arrow  in  illustrates  how  one  passes  from  the 

chosen  “North”  to  the  chosen  “South”  pole  at  the  center  of  each  3-ball. 


Partitioning  ‘Energy’  Between  the  Two  Bath  Modes  Recall  that  (for  3  DoF)  there  are  2  bath  modes. 
For  a  given  trajectory  in  the  NHIM,  all  of  the  ‘energy’  could  be  in  the  first  bath  mode  (resulting  in  a  periodic 
orbit)  or  in  the  second  bath  mode  (resulting  in  another  periodic  orbit),  or  the  energy  could  be  distributed, 
or  partitioned,  between  the  two  bath  modes. 

In  Fig.  36(a)  we  show  the  two  open  balls  that  comprise  our  model  of  as  described  above.  We  choose 
a  trajectory  having  the  property  that  all  of  the  energy  is  in  the  first  bath  mode,  i.e.,  J3  =  0  and  J2  = 
where  is  defined  by  the  energy  equation,  A'nf(0,  0)  =  h  >  0  (see  Sec.  2.2).  This  corresponds  to  a 

periodic  orbit  that  moves  along  the  “vertical”  circle  that  passes  through  the  centers  of  the  two  balls.  As  the 
periodic  orbit  exits  the  North  ball  at  the  top,  it  enters  the  South  ball  from  the  top,  exiting  the  South  ball  at 
the  bottom,  and  entering  the  North  ball  at  the  bottom. 

Suppose  now  that  we  change  this  distribution  of  energy  very  slightly,  by  taking  a  small  amount  of  energy 
out  of  the  first  bath  mode  and  putting  it  into  the  second,  which  we  illustrate  in  Fig.  36(b)  (and  where  we 
leave  it  to  the  reader  to  consider  how  a  trajectory  leaves  one  ball  through  its  surface  and  enters  the  other 
ball).  The  motion  still  stays  close  to  the  original  circle,  but  there  is  now  a  second  circular  (oscillatory)  motion 
“around  the  periodic  orbit”.  Thus  the  motion  now  takes  place  on  a  thin  2-torus.  Any  initial  condition  on 
the  2-torus  corresponds  to  a  quasiperiodic  trajectory  that  remains  on  the  2-torus. 

We  continue  this  procedure  of  taking  energy  from  the  first  bath  mode  and  putting  it  into  the  second  bath 
mode,  until  all  of  the  energy  is  in  the  second  bath  mode  (Fig.  36(e)),  which  results  in  a  periodic  orbit  moving 
in  a  “horizontal”  circle  along  the  surface  of  the  two  balls. 
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Figure  36:  Illustration  of  the  invariant  tori  in  SNfflM(^)j  for  n  =  3  degrees  of  freedom,  (a)  representation 
of  the  periodic  orbit  corresponding  to  local  integrals  of  motion  {I,J2,J'i)  =  (0,  shown  in  both 

hemispheres  of  the  model,  (b)  a  nearby  torus  with  J3  =  e  with  0  <  e  (c)  the  case  J2  =  (d) 

a  torus  with  J3  =  —  e  which  is  close  to  the  periodic  orbit  corresponding  to  (7,^2, -/a)  =  (0,0, 

shown  in  (e). 
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In  this  way  we  see  that  is  foliated  by  a  one-parameter  family  of  2-tori.  Trajectories  on  the  2-tori 
correspond  to  quasiperiodic  motion  with  energy  contained  in  both  bath  modes.  The  “extreme”  cases  corre¬ 
sponding  to  all  the  energy  contained  in  one  of  the  bath  modes  corresponds  to  periodic  orbits  (the  “Lyapunov 
orbits”).  We  illustrate  this  family  of  2-tori  in  Fig.  37. 


Figure  37:  The  family  of  invariant  2-tori  that  fill  the  NHIM,  S^himj  illustrated  according  to  our  method 
for  visualizing  in  three  dimensions.  One  can  think  of  each  torus  as  being  cut  in  half,  with  one  half  in 
the  North  ball  and  the  other  half  in  the  South  ball,  and  the  “edges”  of  the  two  halves  touching  touching 
the  surface  of  the  3-ball,  where  they  are  identified.  In  this  way,  a  quasiperiodic  trajectory  winding  around  a 
2-torus  continually  moves  from  one  3-ball  to  another  as  it  encounters  the  edge  of  the  torus  where  the  torus 
is  cut  into  two  pieces,  which  are  on  the  surfaces  of  the  3-balls  (i.e.  the  “equator”  of  §nhim)' 


5.2  The  NHIM,  Snhim(^)5  forn  >  3  DoF 

The  method  for  visualizing  the  NHIM  for  3  DoF  described  in  the  previous  section  is  not  as  useful  for  n  >  3 
DoF.  Nevertheless,  the  dynamics  on  the  NHIM  have  the  same  characteristics-quasiperiodic  trajectories  with 
the  energy  partitioned  amongst  the  n  —  1  bath  modes.  Hence,  as  mentioned  in  Sec.  4,  for  the  general  DoF 
case,  a  description  of  the  NHIM  can  be  given  in  terms  of  the  geoemtry  of  the  space  of  centre  actions  (or 
“bath  mode”  actions)  J2,  ■  ■  ■ ,  Jn-  This  is  what  we  will  now  describe. 

Classically,  the  space  formed  by  the  centre  actions  is  the  (n  —  l)-dimensional  positive  cone;  equivalently, 
the  positive  orthant  (the  generalization  of  a  quadrant  in  the  plane  or  an  octant  in  three  dimensions  to 
arbitrary  dimensions)  given  by  J2, . . . ,  Jn  >  0  in  In  general,  the  NHIM  is  represented  in  this  space  by 

a  (topological)  (n  —  2)-simplex  having  (n  —  1)  vertices,  one  for  each  centre  action.  The  vertices  are  always 
located  on  the  coordinate  axes;  each  vertex  thus  corresponds  to  (n  —  2)  of  the  integrals  Jk  being  equal  to 
zero.  Note  that  the  projection  of  the  NHIM  to  this  space  of  actions  will  generally  be  a  nonlinear  embedding 
of  a  simplex;  the  embedding  is  only  linear  (that  is,  all  cells  of  the  simplex  are  flat)  in  the  case  where  the 
vector  field  defined  by  Hamilton’s  equations  is  linear.  We  will  denote  this  embedded  simplex  by 

A=[2,3,...,n],  (54) 

using  the  index  of  each  centre  action  to  label  a  vertex.  Within  this  simplex,  each  (fc  —  l)-cell  (for  every 
k  =  1, . . . ,  (n  —  1)),  denoted  [ji, . . . ,  jfe],  represents  a  (k  —  l)-parameter  family  of  invariant  A:-tori  which  we 
denote  by 


.  .  . 
Jlj2---3k 


(55) 


Here  _)i ,,...,  jfe  are  mutually  different  elements  of  {2, . . . ,  n}  which  indicate  the  center  actions  Jj^, . . . ,  Jj^. 
are  different  from  zero.  The  family  is  parameterised,  for  example,  by  the  centre  actions  ,  Jj^-n  with 

the  value  of  Jj^  being  determined  by  the  energy  equation,  as  we  will  describe  below. 
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For  example,  each  vertex  (i.e.,  0-cell),  denoted  [j]  (for  j  =  2,...,n),  represents  a  single  periodic  orbit 
Tj  =  S^,  in  which  all  actions  except  Jj  are  set  to  zero  and,  therefore,  the  remaining  action  is  fixed  at  a  value 
Jj  =  which  solves  the  following  equation  that  corresponds  to  restriction  to  a  single  energy  surface, 

iFNF(/,J2,...,Jn)=h>0,  (56) 

subject  to  setting  Ji  =  0  for  all  £  =  2, ...  ,n,  £  j.  The  periodic  orbit,  Tj,  projects  to  a  single  circle  of 
radius  in  the  (gfc,pfc)-plane  of  the  normal  form  coordinates  and  to  the  origin  in  all  other  normal  form 

planes,  {qi,pe)  for  £^  j. 

Each  edge  [ji,  J2]  (for  ji  ^  j2)  of  the  simplex  represents  a  1-parameter  family  of  invariant  2-tori,  in  which 
all  actions  except  -Jj^  and  have  been  set  to  zero.  Without  loss  of  generality,  we  can  choose  Jj^  as  a 
parameter;  each  value  Jj^  G  [0,  fixes  the  value  of  Jj^  by  virtue  of  the  above  iso-energetic  condition 

subject  to  the  modified  restriction  that 


=  0  for  £  ^  ji,  J2-  (57) 

We  denote  this  1-parameter  family  of  invariant  2-tori  in  by 

By  allowing  k  centre  actions  to  be  non-zero  in  the  NHIM,  corresponding  to  a  (fc  —  l)-cell  [ji, ^2  ■  •  ■  ,jk]  of 
the  simplex,  we  can,  without  loss  of  generality,  take  Jj^,  Jj^, . . . ,  Jj^_^  as  parameters  which  fix  the  value  of 
■Jj^,  via  the  condition 

i7NF(/,J2,...,J„)  =h>0,  (58) 

subject  to  setting 


=  0  for  £  ^  ji , i2,  •  ■  • ,  jk-  (59) 

This  defines  a  {k  —  l)-parameter  family  of  invariant  /c-tori,  denoted 

A  schematic  representation  of  the  projection  of  §NfflM(^)  space  of  centre  actions  J2,  J3, . . . ,  Jn 

for  n-DoF  is  shown  in  Fig.  38  for  (i)  n  =  3,  (ii)  n  =  4,  (iii)  n  =  5.  In  (iii)  we  show  a  schematic  3-dimensional 
projection  of  the  4-dimensional  space  of  centre  actions. 


Directional  flux  across  the  dividing  surface.  In  [Waalkens  &  Wiggins (2004)]  it  is  shown  that  the 
(directional)  flux  across  a  given  iso-energetic  dividing  surface  S^"“^(/i)  (that  is,  the  equal  and  opposite  fluxes 
through  the  forward  and  backward  dividing  surface  hemispheres  s-'^d  may  be  reduced  to 

a  certain  integral  taken  over  §NfflM(^)-  representation  of  the  NHIM  in  in  terms  of  the  centre  integrals 

given  above,  this  flux  corresponds  to  the  (n  —  l)-dimensional  volume  of  the  finite  region  of  the  positive 
orthant  bounded  by  the  (n  —  2)-simplex  representing  the  NHIM  times  (27r)"“^.  That  is  up  to  a  factor,  the 
flux  corresponds  to  the  volume  of  an  (n  —  l)-simplex  formed  by  adding  the  origin  to  the  vertex  set  of  the 
NHIM  in  the  space  of  integrals.  In  Fig.  38  for  (i)  and  (ii),  the  (n  —  l)-dimensional  volume  F  corresponding 
to  the  magnitude  of  the  directional  flux  across  the  dividing  surface  hemispheres  and  i® 

indicated.  Applications  of  this  in  both  the  classical  and  quantum  case  can  be  found  in  [Waalkens  et  al.(2008)]. 

6  Conclusions 

This  paper  is  concerned  with  strategies  and  approaches  for  visualizing  the  phase  space  structures  that  govern 
reaction  dynamics.  In  Section  2  we  reviewed  recent  results  that  develop  the  phase  space  geometry  near 
an  equilibrium  point  of  Hamilton’s  equations  of  saddle  x  centre  x  •  •  •  x  centre  stability  type  (‘saddle’  for 
short)  in  n  degree-of-freedom  Hamiltonian  systems.  The  key  technique  used  was  the  method  of  Poincare- 
Birkhoff  normal  forms.  This  method  provided  a  transformation  (and  its  inverse),  valid  in  a  phase  space 
neighborhood  of  the  equilibrium  point,  to  normal  form  coordinates  in  which  the  phase  space  structures  could 
be  explicitly  realized  (and  then  mapped  back  to  the  original  coordinates  via  the  inverse  of  the  normal  form 
transformation).  The  key  phase  space  structure  was  seen  to  be  a  normally  hyperbolic  invariant  manifold 
(NHIM),  which  is  the  higher  dimensional  analogue  of  a  “saddle  point”.  The  NHIM  is  the  anchor  upon 
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Figure  38:  Schematic  representation  of  the  projection  of  space  of  centre  actions, 

J2,  J3,  ■  ■  • ,  Jn,  for  n  DoF  giving  a  nonlinear  embedding  of  a  simplex  with  (n  —  1)  vertices,  [j],  j  =  2, . . . ,  n, 
which  encloses  a  volume  proportional  to  the  directional  flux,  F. 


which  a  dividing  surface  having  the  “no-recrossing”  property  and  minimal  flux  is  constructed.  Moreover,  the 
stable  and  unstable  manifolds  of  the  NHIM  are  the  “conduits”  for  forward  and  backward  reactive  trajectories 
(trajectories  evolving  from  reactants  to  products,  and  vice-versa).  More  precisely,  one  the  two  branches  of 
the  stable  manifold  can  be  glued  together  with  one  of  the  two  branches  of  the  unstable  manifold  to  form  a 
cylinder  (denoted  the  forward  reactive  cylinder  in  this  paper)  that  enclose  all  forward  reactive  trajectories; 
and  the  other  branches  of  the  stable  and  unstable  manifolds  can  be  glued  together  to  form  a  cylinder  (the 
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backward  reactive  cylinder)  which  encloses  all  backward  reactive  trajectories.  In  the  case  where  a  generic 
non-resonance  condition  holds  on  the  pure  imaginary  eigenvalues  associated  with  the  center  directions  it 
follows  that  in  the  normal  form  coordinates  the  Hamiltonian  (when  the  series  representing  the  normal  form 
Hamiltonian  is  truncated  at  some  finite  order)  has  n  integrals  of  the  motion  and  thus  is  integrable.  This 
integrable  structure  allows  for  a  natural  decoupling  of  the  dynamics  into  a  single  reactive  mode  and  n  —  1  bath 
modes.  Moreover,  the  integrable  structure  provides  a  foliation  of  the  NHIM  by  a  family  of  lower  dimensional 
tori  and  a  foliation  of  a  neighborhood  of  the  equilibrium  point  into  Lagrangian  cylinders. 

The  “decoupling”  of  the  vector  field  and  the  associated  dynamics  masks  the  complete  picture  of  reaction 
governed  by  the  phase  space  structures  in  the  energy  surface.  In  Section  3  we  provide  three  different  2 
degree-of-freedom  models  that  allow  us  to  visualize  the  energy  surface,  the  NHIM  and  its  stable  and  unstable 
manifolds,  the  dividing  surface,  and  the  Lagrangian  cylinders.  These  give  a  complete  picture  of  how  the 
phase  space  structures  govern  the  dynamics  of  trajectories  near  a  saddle.  By  contrast,  the  role  of  the  phase 
space  structures  is  typically  obscured  in  their  common  representation  as  projections  to  conhguration.  For 
example,  the  projection  of  the  volume  enclosed  by  the  forward  and  backward  cylinders  might  not  lie  inside  of 
the  projection  of  the  forward  and  backward  cylinders,  and  as  a  consequence  forward  and  backward  reactive 
trajectories  might  seemingly  leave  the  forward  and  backward  reactive  cylinders  (for  an  example,  see  the  study 
of  the  planar  Hill  problem  in  [Waalkens  et  al.  (2005b)]). 

For  the  general  case  of  n  >  2  degrees  of  freedom  systems,  we  use  the  integrals  alone  to  provide  a  model 
for  visualizing  the  phase  space  structures  (see  Section  4  ).  This  model  does  not  contain  explicit  information 
on  trajectories,  but  it  can  be  easily  related  to  the  behavior  of  trajectories  by  utilizing  the  foliation  of  the 
neighborhood  of  the  NHIM  by  Lagrangian  cylinders. 

We  moreover  discussed  in  some  detail  the  foliation  of  the  NHIM  for  3  degrees  of  freedom.  Here  the  NHIM 
is  a  threedimensional  sphere,  that  is  foliated  by  a  one-parameter  family  of  invariant  2-tori.  We  visualized 
this  foliation  using  the  3-ball  model  of  (see  Section  5).  For  the  general  case  of  n  >  3  degrees  of  freedom, 
we  presented  the  NHIM  as  a  simplex  in  the  space  of  the  integrals.  This  is  of  particular  importance,  since, 
up  to  a  prefactor,  the  volume  enclosed  by  the  simplex  in  the  positive  orthant  of  the  integral  space  gives  the 
directional  flux. 

We  conclude  by  emphasizing  that  we  focused  the  discussion  of  the  phase  space  structures  and  their 
visualization  to  the  neighborhood  of  the  saddle.  In  fact,  depending  on  the  global  topology  of  the  energy 
surface,  the  local  partitioning  of  the  energy  surface  by  the  dividing  surface  constructed  from  the  NHIM  into 
a  reactants  and  a  products  component  might  not  extend  to  a  global  partitioning.  For  example,  there  could 
be  more  than  one  saddle  point  which  controls  the  access  to  and  escape  from  a  certain  phase  space  region 
which  requires  a  more  careful  dehnition  of  a  reactants  and  a  products  region.  Moreover,  in  this  paper  we 
did  not  discuss  the  global  geometry  of  the  stable  and  unstable  manifolds  of  the  NHIM.  In  fact  the  stable 
and  unstable  manifolds  might  extend  to  regions  far  away  from  the  neighborhood  of  the  saddle.  Even  if  there 
is  only  a  single  saddle  point,  the  global  structure  might  be  very  complicated  since  the  stable  and  unstable 
manifolds  might  wrap  around  in  a  complicated  manner  and  intersect  each  other  to  form  a  highdimensional 
homoclinic  tangle.  If  there  are  several  saddles  the  stable  and  unstable  manifolds  of  NHIMs  of  different  saddles 
might  intersect.  Since  the  stable  and  unstable  manifolds  are  of  codimension  1  in  the  energy  surface  they  still 
act  as  impenetrable  barriers  globally,  and  also  control  the  global  dynamics  in  a  possibly  complicated  way. 
The  details  of  the  dynamics  then  strongly  depends  on  the  system  under  consideration.  For  a  few  examples, 
we  refer  to  [Waalkens  et  al. (2004a),  Waalkens  et  al. (2004b),  Waalkens  et  al. (2005a),  Waalkens  et  al. (2005b), 
Waalkens  et  al.  (2005c)]. 
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